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Abstract of the Dissertation 

Study of Open Charm Production in p+p 
CoUisions at = 200 GeV 

by 

Sergey Butsyk 
Doctor of Philosophy 
in 
Physics 
Stony Brook University 
2005 



The PHENIX experiment at the Relativistic Heavy Ion Colhder 
(RHIC) with its unique electron identification system enables us 
to perform high precision measurements of electron yields. By 
measuring electron production at high transverse momentum, we 
can disentangle the contribution of electrons originating from semi- 
leptonic decays of heavy quarks (charm or bottom) from the less 
interesting "photonic" decay modes of light mesons. D/B mesons 
carry single heavy valence quarks and are usually referred to as 
"Open Charm" and "Open Bottom" particles, differentiating them 
from Closed Flavor particles such as J/ip and Y mesons. Due to 
the large mass of the heavy quarks, their production mechanisms 
can be adequately explained by perturbative QCD (pQCD) theory. 

This dissertation presents the measurement of electrons from 
heavy fiavor decays in proton + proton collisions at RHIC at colli- 
sion energy ^/s = 200 GeV over a wide range of transverse moment 
(0.4 < pt < 5 GeV/c). Two independent analysis techniques of 
signal extraction were performed. The "Cocktail" subtraction is 
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based on the calculation and subtraction of the expected "photon 
- related" electron background based upon measured yields of light 
mesons. The "Converter" subtraction is based upon a direct mea- 
surement of photon yields achieved introducing additional material 
in the PHENIX acceptance and deducing the photon abundance by 
measuring the increase in electron yield. This is the first measure- 
ment of the Open Charm crossection at this collision energy and it 
is an important baseline measurement for comparison with nucleus 
-I- nucleus coUisions. The modification of Open Charm production 
in heavy ion collisions compared to the presented p + p result can 
be used to study the final state interaction of the heavy quarks 
with hot dense matter inside the collisions. The results of the 
Open Charm measurements are compared to current pQCD pre- 
dictions both in Leading Order (LO) 0{al) and Next-to-Leading 
Order (NLO) 0{al). The final result for the total Open Charm 
crossection is dec = (0.920 ± 0.148(stat) it 0.524(sys)) mb, in good 
agreement with the published STAR experiment measurements of 
Open Charm crossection in d + Au colliding system. 
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Chapter 1 
Introduction 



Heavy quarks are a unique tool in High Energy and Nuclear Physics. Dis- 
covered in 1974, the charm quark has attracted the attention of both exper- 
imentalists and theoreticians. The production of particles with Open Charm 
(carrying a single c or c quark, such as D-mesons and baryons) and Closed 
Charm (cc mesons) is a rare process purely from energy considerations, be- 
cause the charm mass, mc ~ 1.2 — 1.5 GeV, is so high. This energy scale allows 
one to use (with some caution, of course) perturbative Quantum Chromody- 
namics (pQCD) - a theoretical model of the strong color field interaction to 
describe the production mechanisms and rates of Open Charm. 

Open Charm particles are produced through the fragmentation of cc pairs. 
Creation of those pairs is intrinsically sensitive to the initial state of the par- 
tonic system. pQCD predicts that the primary production mechanism in 
Leading Order approximation is "gluon fusion". Thus, the production rate 
of charm pairs is proportional to the initial gluon density inside the collision. 
Open Charm production in in heavy ion collisions at the Relativistic Heavy 
Ion Collider (RHIC) can thereby lead us to understanding how the initial prop- 
erties of the matter change with the size of the coUiding system. RHIC raised 
the High Energy Nuclear physics to the next level, enabling us to study with 
the same detector setup the full variety of collisions - p + p, d + Au, Au + Au 
at \/i = 200 GeV. Detailed comparison of Open Charm production in those 
collisions will help us to answer the fundamental questions: 

• What are the mechanisms governing charm production? 

The p-\- p collision is a reference measurement, which should be exactly 
explainable by perturbative QCD. Production of Open Charm at lower 
energies has always been a serious confirmation of pQCD's applicability 
for heavy quarks. So far, exact calculations from pQCD exist to second 
order (NLO) [2, 3, 4, 5, 6] and describe the lower energy data within 
the theoretical uncertainties. It may well be that at our high energy, the 
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theory cannot neglect these higher order partonic interaction processes 
in order to describe the data. The theory of hadronization of the c and c 
is presently at a rudimentary stage and in order to constrain the differ- 
ent theoretical predictions, we need to study p + p collisions at different 
rapidity ranges. 

• How does the production sensitive to additional "cold" matter in the 
collision? 

A comparison of p + p production with d + Au will measure the mod- 
ification of the gluon structure function in "cold" hadronic matter. In 
this case, the ability to perform rapidity scan is essential to study the 
X dependence of the gluon modification factor that was clearly observed 
in case of light quark mesons. 

• What happens with the charm quark production in "hot" dense matter 
of Au + Au collisions? 

Direct comparison of charm production crossection in Au -\- Au with 
d -\- Au and p + p will answer this question. So far there are a lot of dif- 
ferent predictions, expecting both enhancement due to the initial state 
thermal production and suppression due to the media induced gluon ra- 
diation. So far, we observed a significant suppression for the light quarks, 
in case of the heavy quark this effect may be significantly reduced due to 
the suppression of gluon radiation at small angles for the charm quark. 

There is a lot of very interesting new physics that can be discovered 
if this programm is accomplished. In this thesis we begin the program by 
accomplishing the baseline measurement of Open Charm production in p -|- p 
collisions at mid rapidity. 

There are two experiments at RHIC capable of performing this task, each 
with certain advantages and disadvantages - PHENIX and STAR. 

PHENIX was a pioneer in Open Charm measurements at RHIC. Built 
with a strong accent towards the electron and muon identification, PHENIX 
successfully measured Heavy Flavor production at ^/s — 130 GeV [50]. Fig. 1.1 
shows the invariant multiplicity for the electrons, originated from the semi- 
leptonic decay of the charm and bottom quarks c{h) —>■ eX. 

STAR also has elaborated electron identification using Electromagnetic 
Calorimeter, Time-of-Flight Detector and Time Projection Chamber. They 
were able to perform similar semi-leptonic measurement in p + p and d + Au 
collisions at ^/s = 200 GeV [72]. But STAR has also been able to directly 
measure decays of D mesons into hadronic channels. This is a more direct 
measurement of the Open Charm production then measuring decay leptons. 
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Figure 1.1: Invariant crossection for 
"Non-photonic" electrons at ^/s = 130 
GeV at PHENIX [50]. 



Figure 1.2: Mass spectrum of Ktt 
pairs in d+ Au collisions with clear D° 
peak [72] 



Fig. 1.2 shows the reconstruction of Kn decay channel by STAR col- 

laboration. It is doubtful that these measurements can be also done in case of 
Au + Au collisions. Direct measurement of hadronic decays of Open Charm 
mesons in PHENIX is very hard due to the limited acceptance for these decays. 

My dissertation sets the cornerstone to the foundation of the Heavy Fla- 
vor measurements at energy ^/s = 200 GeV - I will present the experimental 
measurement of Open Charm production at mid-rapidity in PHENIX at RHIC 
from single electron channel (c — eX). The paper is organized in the following 
way: Chapter 2 presents the current theoretical aspects of the heavy quark 
production and fragmentation. The experimental setup, used in the analysis 
is described in Chapter 3. Chapter 4 is dedicated to data reduction towards 
the final Open Charm electron crossection. A discussion and theoretical com- 
parison is summarized in Chapter 5 and final conclusions and a future outlook 
are presented in Chapter 6. 
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Chapter 2 



Theory 



2.1 Quantum Chromo dynamics 



Quantum ChromoDynamics (QCD) is the set of physical laws, describing 
our current understanding of strong interactions between the quarks and glu- 
ons. QCD was developed as an extension of Quantum Electrodynamics via the 
imposition of a local SU (3) symmetry of rotation in color space. The Standard 
Model couples the SU (3) strong interaction with the unified electro- week inter- 
action making the overall symmetry group SU (3) x SU (2) x f/(l)y. The Stan- 
dard Model describes the large variety of mesonic and baryonic states by the 
existence of the deeper level of elementary constituents of matter: quarks [8]. 
All the variety of particles that can be found in Particle Data Tables [9] can 
be explained in terms of six spin-| quark flavors: 



Q = H 


u c t 


Q = -k 


d s b 



The meson state in this case consists of quark-antiquark (gg) and baryonic 
state is a combination of three quarks {qqq) or antiquarks (ggg), this model 
gives us a unique correspondence between the observed particle and its state. 

However, the early quark picture faced a problems with Fermi-Dirac 
statistics, for example A++ baryon (J = |) should have been explained as 
a combination of m | m | m | quarks which contradicts Fermi statistics of 
spin- 1 particles. This problem was solved by introducing a new quantum 
number [10], color, such that each species of quark can have Nc = 3 differ- 
ent colors. In this picture mesons and baryons are described as color-singlet 
combinations of colored quarks. 



B 





Qaqp) 



(2.1) 




p 



Figure 2.1: Deep Inelastic e p — > e X scattering. [1] 



States, with color content not equal to zero were never observed, which led 
the theorists to postulate the hypothesis that all states are color-singlets. This 
statement is known as "confinement hypotheses" implying directly the 
non-existence of free quark colored states. 

The next valuable piece to the construction of the QCD theory was ob- 
tained from "deep inelastic scattering" (DIS) of high energy electrons on 
the proton. Fig. 2.1 shows the interaction of of the electron with the content 
of the proton. 

At higher energies, the virtual photon is sensitive to smaller distances 
and we are thereby able to resolve the structure of the proton by increasing 

the energy of the electron beam. Denoting the photon momentum transfer 
as = 4:EE'sin'^{9/2), where E and E' are incident and scattered electron 
energies. 

From the studies of scattering crossections, it was found that at large 
values of a sizable continuum contribution to the crossection still persists, 
suggesting an existence of pointlike objects inside the proton. Those point- 
like spin- 1 constituents were called partons. It was also found that proton 
structure functions only depends on the transferred momentum and the 
relative contribution of parton momentum to the total proton momentum 
X [11]. This surprising discovery of pointlike partons inside the proton, ini- 
tially contradicting previously explained confinement of the quarks by a strong 
color force leads to an fundamental property of strong interactions known as 
" asymptotic freedom" - the interaction between quarks should become weaker 
at short distances, so that at high Q"^ transfer, quarks behave as free particles. 
The interaction in qq system grows as we try to separate the quarks. At some 
point the energy of the string become more then 2mQ, where itlq is the mass 
of another quark. It becomes energetically favorable to split the string into 
two qQ and Qq mesons. By increasing the energy, we can create more and 
more colorless meson states. Otherwise, if we try to approach two quarks, the 
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potential loses strength and the qnarks start behave like free particles. 

The field carriers in QED theory are photons, and the QED Lagrangian 
describes the interaction of a Dirac fermion with the electromagnetic field 
A^. By analogy, the QCD Lagrangian can be constructed. This Lagrangian 
will describe the interaction of g^g^ quarks of color a and flavor /. The free 
Lagrangian is given by Eq. 2.2 [1]. 

Lo^J2^f (^^''^M - mf)qf (2.2) 
/ 

Propagation of the quark wave function requires the Lagrangian to be 
invariant under local SU{3)c color space transformation. This invariance in- 
troduce 3^ — 1 = 8 component gauge bosons G'^{x), those particles are called 
gluons. The total invariant Lagrangian of QCD can be written as follows [1]: 



f 

f V ^ / q/3 

_ I fabc _ Q.^^^ ^pc _ 9l ja,cj^^ G^.G^^^Gpl 

The first line presents the kinetic energy of quark and gluon field. The 
second line describes the color interaction between gluon and quark, where A" 
are color field SU{3)c matrices. The last line shows gluon self-interactions, 
exhibiting non-abehan properties of the color field^. The coupling constant 
Qs is the universal strong coupling constant. Examples of qq gg scattering 
Feynman Diagrams allowed in QCD are shown in Fig. 2.2 [1]. 



q G 




q G 

a) b) 0) 



Figure 2.2: Tree-level Feynman diagrams contributing to qq — > gg. [1] 



Gluon field self-interaction is a striking difference of QCD, compared to QED, 
where photon can not spUt into two photons 
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2.2 Renormalization of QCD 



Renormalization techniques have been widely developed in QED as a use- 
ful mechanism to absorb the loop diagrams that cause logarithmic divergences. 
Ultraviolet logarithmic divergences appear in QED if the photon polarizes 
the vacuum, creating a virtual e"'"e~pair. The calculated crossection for this 
loop process produces an infinity which is "non-physical". In order to solve 
this problem, we assume (see Fig. 2.3) that for e^e^intcractions processes, 
the effective coupling constant that we measure is changing due to vacuum 
polarization. The "bare" coupling constant in QED has been known since 
Thompson experiments as ao — e^/ (An) ^ 1/137. Assuming that the couphng 
may vary as a function of photon transferred energy due to appearance of 
the loops leads to very powerful conclusion that the we can actually measure 
renormalized running coupling a^^^{Q'^) which can be written as: 



where /Sq is a solution of differential equation (2.5) and for one loop con- 
tribution is equal to = |. 



Qq is an arbitrary scale, at which we can neglect the contributions from 
the radiative corrections. We can take it as Ql ^ which produces classical 
value for the coupling constant a{Ql) — ao — 1/137. 

The fact that f3o > leads to the conclusion that alpha is growing with 
the increase of the energy transferred and, i.e. the effective charge of the elec- 
tron is decreasing with the distances between the charges. This is intuitively 
understandable as the produced virtual e"'"e~pairs acts like a screen around 
the charge in the polarized dielectric medium. 




(2.4) 





(2.5) 



e" 



e" 




e- 



Figure 2.3: Photon self-energy contribution to e e scattering. [1] 



Figure 2.4: Feynman diagrams contributing to the renormalization of the 
strong couphng. [1] 



A similar renormahzation approach can be apphed to the QCD theory, 
but in the case of QCD we need to take into account gluon self-interaction 
diagrams. Processes, contributing to the renormalization of the scale are shown 
in Fig. 2.4. The final equation for the QCD running coupling gives identical 
result to QED Eq. 2.4, but the /3 function will be drastically different [12]. 



A ^ '-^^i^ (2.6) 

where Nf is a number of flavors, Nc is the number of colors. The first term 
is due to the creation of qq loops of different fiavor and is similar to the QED 
value. The second, negative contribution in Eq. 2.6 corresponds to gluonic 
self-interaction loops and is proportional to number of colors. The existence 
of the second term makes the factor Pq negative for Nf < 16 this leads to the 
important result that that at short distances (large Q^) as{Q^) — > 0. This is 
direct confirmation of asymptotic freedom . 

In the one loop approximation we can write the strong interaction running 
coupling as (solving Eq. 2.5 for a): 



where Aqcd comes as an integration parameter at an energy scale at which 
the QCD coupling constant starts to exhibit non-perturbative properties. If 
the momentum transfer Q is much larger then Aqcd, we are in asymptotic 
freedom regime and as{Q^) — > allows us to use the perturbative QCD 
(pQCD) technique. In the case of energy transfer, comparable to Aqcd, Eq.2.7 
blows-up the strong coupling constant making all orders diagrams to be of 
similar contribution and the pQCD approach cannot be used. 
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It is necessary to add that there are several theoretical approaches to 
pQCD renormalization (referred to as the renormalization scheme). Currently, 
the most adopted and cited approach is use of MS - modified minimal sub- 
traction scheme [13]. 

In MS case the next orders of two loop corrections we have following 
expression for next term of (3 function deconvolution Eq. 2.5 [14]. 

51 1 Q 

A = -X + li^/ (2.8) 

For Nf < 8, /3i < which even stronger justifies the use of asymptotic 
freedom. Fig. 2.5 shows the dependence of on the energy scale Q for the 
experimental data comparison with pQCD theory [15]. It is also clearly seen 
from this figure that Aqcd is on the order of ~ 220 MeV and should be 
calculated for the expected number of flavors (Ag^^ assumes Nf = 5). 
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Figure 2.5: Measurements of as a function of energy scale [15]. 



2.3 Heavy quark production in QCD 

In this section I discuss the production mechanisms of heavy quark par- 
ticles in p + p coUisions using the perturbative QCD formahsm, described in 
previous Section. 

The main advantage of the heavy quark measurements is that the initial 
content of the heavy flavor in the proton is negligibly small (there are some 
papers though that assume the intrinsic heavy quark content of the nucleon to 
be non-zero [16, 17] but we will not include this aspect into our considerations). 
Then, the heavy quarks should be produced in hard partonic scattering. The 
big advantage of the heavy quark production is that the energy scale for the 
production of charm quarks ~ is significantly higher then Aqcd- This 
gives us a coupling constant of the order of as ~ (0.3 — 0.5) (see Fig. 2.4), 
which is small enough to utilize perturbative theory. pQCD should give even 
more accurate results for the higher mass bottom quark. 

2.4 Hard scattering processes 

In the partonic model we can assume that the colliding protons consist 
of a collection of partons (quarks, antiquarks, gluons) that interact with the 
content of the other proton. In Leading Order pQCD, hard scattering of the 
partons can create heavy quark- antiquark pairs QQ^ by the diagrams, shown in 
Fig. 2.6a,b). LO sub-processes are usually called "gluon fusion" {gg — > QQ) 
and " quEirk-antiquark annihilation" (qq — > QQ). 

The general perturbative extension for the total crossection of quark pair 
production on the partonic level can be expressed by the following equa- 
tion [19]. 

-rAs.ml.,l) = {^--s{,l)TY. /r(^) ^ (2-9) 

k=o 1=0 y^Qj 

where s = X1X2S is the partonic energy in center of mass. The useful di- 
mensionless parameter r] = s/AttIq — 1 [20] reflects the threshold production of 
the heavy quark (v^ should be at least 2rnQ to create a quark-antiquark pair). 
i and j are the partonic indexes corresponding to each particular Feynman di- 
agram on the Fig. 2. 6. fij'^\r]) is a dimensionless scaling function representing 



index refers to heavy quarks 




Figure 2.6: LO and most important NLO heavy quark production diagrams. 
LO - a) "gluon fusion" b) "quark-antiquark annihilation" NLO - c) Pair cre- 
ation with gluon emission in output channel d) "flavor excitation" e) "gluon 
splitting" f) "gluon splitting but of "flavor excitation" character [18] 

the amplitude of a given partonic scattering diagram. Index k shows the order 
of the process diagram, k = corresponds to LO, one can see that in the 
leading order the total crossection is a function of a^. Next order of pertur- 
bation k = 1 corresponds to Next-Leading-Order (NLO) contributions (some 
of the NLO diagrams are shown in Fig. 2.6c-d)). The NLO crossection is of 
the order of and in principle are of the next order of smallness, unfortu- 
nately, the NLO crossection develops a logarithms term for I = 1 which can 
ln{-^) be make NLO crossection of the order and even higher then the LO 

contribution. Case of /c = 2 is usually referred to as Next-to-Next Leading 
Order (NNLO) so far the exact crossection for NNLO heavy quark production 
has not been calculated because of a large number of contributing diagrams 
and mathematical complications in resummation procedure. There have been 
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theoretical attempts to calculate the NNLO crossection contribution near pro- 
duction threshold [20]. 

in Eq. 2.9 is called renormalization scale and usually assumed to 
be proportional to or mg ^Vtq depending on the model. The propor- 
tionality coefficient sc^ in = [scu • mqf is a variable parameter in QCD 
and usually assumed to be in the range of scr G (|; 2). Low values of scj? = \ 
for charm quark creation produces scales on the order of < 1 GeV^ and 
from Fig. 2.5 one can see that a.s becomes quite large to utilize perturbative 
technique. 

The next general step in QCD is to assume that we know the distribu- 
tion of the partons inside of the proton. This distribution is called peirton 
distribution function (PDF) [21, 22] ff(x,ijLp) described in terms of Feyn- 
man variable x and the momentum transfer scale jJ/jp- X IS Si relative fraction 
of partonic momentum to the total momentum of the hadron Xi = . 

Pz max 

Scale dependence of the PDF is described by DGLAP equations [23, 24, 25] 
and the shape of the parton distribution is derived from comparison of pQCD 
predictions to experimental measurements. 

Fig. 2.7 shows the most striking agreement of pQCD prediction and the 
results of 7r° and direct-7 production in p+Be coUisions from £'706 experiment. 
The results show the existence of non-zero smearing of the transverse (/ct) ~ 1 
GeV of the parton (so called intrinsic kx). Experimental results helped to 
constrain the gluon distribution function with very good accuracy. There is 
a big variety of PDF functions developed by different theory groups. At the 
current moment, the CTEQ group's PDF set [22] seems to be the most accurate 
in describing existing experimental data for the structure of the proton. 

Taking the assumption of the parton distribution, we can write the total 
crossection for Heavy Flavor production [19] 



'pp\ 



(2.10) 



where aij{TS,mQ, /i'^) is the partial partonic crossection (Eq. 2.9), Hp is 
a momentum transfer scale (factorization scale) of the PDF factorization. 
The usual assumption for pQCD calculations is to expect the factorization and 
renormalization scales to be similar fij^ = /j,^ = jj,^ since all PDFs uses this 
assumption a priory in calculations [19]. 
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Figure 2.7: Comparison of the photon and direct-7 measured by E706 ex- 
periment in p + Be colhsions at ^/s = 38.8 GeV with pQCD predictions for 
different values of (/ct) smearing (left panel). Comparison of gluon structure 
function for most resent PDF sets (right panel). [26] 



2.5 QQ fragmentation 

After the quark- ant iquark pair is created, each of the produced heavy 
quarks picks up a light quark from the proton in order to create a color singlet 
object (string). The string is fragmented into hadrons using a phenomenolog- 
ically determined fragmentation function D{z) where z is the momentum 
fraction of the quark, carried by the hadron. D{z) determines the probability 
of producing hadron with given momentum. 

By default, PYTHIA (a widely used pQCD Monte Carlo program) uses 
the Lund fragmentation function [27] modified by Bowler [28]: 



(1 — zY 
D{z)lund oc exp 



z 



(2.11) 



where = m\ + p\ is the transverse mass of the hadron. The default 
parameter values in PYTHIA are a = 0.3 and b = 0.58 GeV~^. 
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A different parametrization of the fragmentation function was proposed 
by Peterson [29] 

D{zU^.rs^ oc (2-12) 

The central value for e in D-meson decay usually assumes to be e = 
0.06 [30, 31] in case of LO theory fits to the data, this parameter is much 
smaller if NLO fits are used for parameter extraction (e = 0.01 — 0.02) [32]. 

Now we have everything to calculate the single differential crossection for 
the produced Heavy Flavor hadrons. The non-perturbative hadron production 
can be obtained using the factorization theorem: 

^ = J^PTd'^ ^(^) S{PT - zp^) (2.13) 

where is single differential crossection for heavy quark. 

Heavy Flavor hadroproduction is currently a subject to significant uncer- 
tainty. The most accurate approach to deriving the fragmentation function 
is to use the Mellin transforms of the D{z) and obtain the momenta of this 
transform from the experimental data [33, 34, 35]. Mellin transformation is a 
decomposition of the function on the powers of N according to the following 
formula [34]: 



Dn = 



= J D{z)z''^ (2.14) 



Assuming the quark px distributions have power law behavior ^ = ^ 
in the neighborhood of some pt, one can immediately find from Eq. 2.13 that 

dfjH f A A 

-— = dpTdz — D{z) 5{pT - zpr) = —D^ (2.15) 
dpT J Pt Pt 

Thus, the hadronic crossection is given by the product of n*'* moment 
of the fragmentation function. Recent results from the ALEPH collaboration 
on the bottom meson energy transfer indicates that Peterson fragmentation 
form does not describe the calculated momenta of the fragmentation function 
(see Fig. 2.8). More studies of fragmentation function form for Heavy Flavor 
mesons need to be performed to remove this apparent uncertainty. 
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Figure 2.8: Momenta of measured B meson fragmentation function compared 
with pQCD NLL calculations with different assumptions for D{z). The solid 
line is one-parameter fit to the second momentum [34]. 

2.6 Heavy Flavor meson decay 

Heavy fiavor semi-leptonic decay diagrams for D and B mesons are shown 
in Fig. 2.9. This weak decay channel should be calculable using the standard 
technique, developed for decay of free mesons {spectator model) [36]. 

The total semi-leptonic decay rate can be written as 

" S / ^^^^^(^ + y-Xm) e{x^ -x-y + xy)^ $^|M«|2 (2.16) 

where x = 2Ee/mQ and y = 2Ey/mQ are the rescaled energies of charged 
and neutral lepton in the heavy quark rest frame. Xm = 1 — i^q/^Q)"^ is 
kinematic limit to the energy transfer. 

Matrix decay elements equal to: 

^\M^\^ = 64Gl\V,s\^ c ■ e+ s ■ u ^IM^ = 64G||KfeP & ■ c ■ e+ (2.17) 

where c, b, s, e^, v are corresponding 4- momenta of the decay particles. 
|V^6| are elements of Cabibbo-Kobayashi-Maskawa (CKM) flavor mixing 
matrix. Gp = 1.16637 • 10~^ GeV~^ is Fermi constant of weak interaction. 
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From Eq. 2.16 we can now calculate the differential semi-leptonic decay 
rate for the Heavy Flavor meson. 




Figure 2.9: Semi-leptonic decay diagrams for Open Charm and Open Bottom 
meson semi-leptonic decay. 



dr 



dx 



dx 

where To{mQ) 



\VcbfTo{mb) 



'12x{xm - xf 

(1-X)3 



(2.18) 



{6 — 6x + X Xm + — 3Xm) 



q2 ,^^^5 

-^§2^- -^^S- 2-10 shows the resulting decay rate as a function 
of relative electron energy for c — > s, c — > d, 6 — > c, 6 — > it semi-leptonic 

decays. 

Thus, pQCD theory provides us with the solid prescription to calculate 
the Heavy Flavor production, hadronization and decay. 




Figure 2.10: Semi-leptonic decay rate as a function of 2E^/Mq for Open Charm 
(left panel) and Open Bottom (right panel) [36]. 



Chapter 3 

PHENIX Experiment 



3.1 PHENIX Detector 

PHENIX is one of the most advanced detector systems in High Energy 
Nuclear Physics. It consists of 14 independent detector subsystems collect- 
ing the information about the particles produced in Heavy Ion collisions. 
The PHENIX tracking system provides accurate measurements of charged 
particle tracks and particle identification (electron, muon, hadron identifica- 
tion) over wide range of momentum. Together with charged particle tracking, 
PHENIX also performs a high precision measurements of photons via a large 
area Electro-Magnetic Calorimeter. A sophisticated Level-1 and Level-2 trig- 
ger system, enables a unique ability to address specific rare physics events 
recording at high luminosity {J/'ip leptonic decays, for example, by trigger- 
ing on a pair of highly energetic electron candidates) . All those factors make 
PHENIX both the most challenging and the most capable detector in RHIC 
to study Heavy Flavor particle productions. 

The detectors in PHENIX are grouped into three major categories by 
design and specific physics tasks they are intended to perform: 

• Global Trigger Detectors 

• Central Arm Detectors 

• Muon Arm Detectors 

The PHENIX detector if shown in Figure 3.1. It is symmetric around mid- 
rapidity with the interaction point positioned in the center of the magnets, 
commonly referred to as the Interaction Region (IR). Two Central Arm detec- 
tors (referred to as East and West) are located around the interaction point and 
cover the rapidity range —0.35 < rj < 0.35. Each Central Arm covers 90° in 
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Table 3.1: Summary of the PHENIX Detector Subsystems [37]. 



Element 


Ar/ 


A(f) 


Purpose and Special Features 


Magnet: central (CM) 


±0.35 


360° 


Up to 1.15 T-m. 


muon (MMS) 


-1.1 to -2.2 


360° 


0.72 T-m for 77 = 2 


muon (MMN) 


1.1 to 2.4 


360° 


0.72 T-m for = 2 


Silicon (MVD) 


±2.6 


360° 


d'^N/dr]d(j), precise vertex, 








reaction plane dctcrmmation 


Beam-beam (BBC) 


±(3.1 to 3.9) 


360° 


Start timing, fast vertex. 


NTC 


±(1 to 2) 


320° 


Extend coverage of BBC for 








1 A 

p-p and p-A. 


ZDC 


±2 mrad 


360° 


Minimum bias trigger. 


Driit chambers (UCj 


±0.35 


90° X 2 


1 J 1 

Good momentum and mass 








resolution. 








Am/m = 0.4% at m = IGeV. 


Pad chambers (PC) 


±0.35 


90° X 2 


Pattern recognition, tracking 








C 1 11* J * 

tor nonbend direction. 


TEC 


±0.35 


90° 


Pattern recognition, ah /ax. 


RICH 


±0.35 


90° X 2 


Electron identification. 


ToF 


±0.35 


45° 


Good hadron identification. 








a <100 ps. 


TO 


±0.35 


45° 


Improve ToF timing for p-p 








and p-A. 


PbSc EMCal 


±0.35 


90°+45° 


For both calorimeters, photon 








and electron detection. 


PbCrl EMCal 


±0.35 


45 


Good e I'K separation at 








•p > \ GeV/c by EM shower 








and p < 0.35 GeV/c by ToF. 








/ti^ separation up to 1 








GeV/c by ToF. 


fjL tracker: (//TS) 


-1.15 to -2.25 


360° 


Tracking for muons. 


(/xTN) 


1.15 to 2.44 


360° 


Muon tracker north installed 








for year-3 


/X identifier: (/xIDS) 


-1.15 to -2.25 


360° 


Steel absorbers and larocci 








tubes for 


(//IDN) 


1.15 to 2.44 


360° 


muon/hadron separation. 
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and populated with the particle detectors for precise measurement of electrons, 
photons, and charged hadrons. Two Muon Arm detector subsystems (referred 
to as South and North) measure the muon and decay hadrons yield at for- 
ward and backward rapidity region [t] G (—2.25, —1.15) and rj G (1.15,2.44)). 
Global Trigger detectors provide the Level-1 trigger information for the colli- 
sion by measuring the time and position of the interaction vertex with high 
precision. 

Detector parameters and performance, rapidity and azimuthal angle (0) 
coverage, of each subsystem are summarized in Table 3.1. 



PHENIX Detector - Second Year Physics Run 




West Beam View 




Figure 3.1: Layout for PHENIX Experiment. 



Wc define a coordinate system relative to the beam axis. The origin is 
located in the center of the IR with the Z-axis pointing along the beam line 
from South to North. The X-axis is points horizontally from East to West and 
the y-axis points upwards, making a proper right-handed coordinate system. 

3.2 Global Detectors 

Nuclear collisions are characterized by their impact parameter vector 
(both magnitude and direction) . Both of these parameters can be determined 
with reasonable collision on a event-by-event basis. Nearly all physics mea- 
surements are then studied to determine their variation with respect to these 
so-called "global" characteristics of the event. 

"Global" detectors measure global characteristics of the collisions and is- 
sue online trigger decision to read-out the information. In PHENIX this task 
has been principally accomplished by two subsystems: 
the Beam — Beam Counters (BBC) and the Zero Degree Calorimeters 
(ZDC). The BBC consists of a set of two (South and North) fast trigger coun- 
ters providing information about Z vertex position and collision time with 
respect to the RHIC beam crossing clock {to bbc)- The ZDCs are a set of 
two Tungsten calorimeters (South and North) located far from the collision 
point and providing additional information about the impact parameter of the 
collision. The third global detector is the Multiplicity Vertex Detector 
(MVD), a silicon barrel detector designed to measure (P N / drj d(j) distribution 
for charged particles near the collision point. This detector is not yet func- 
tioning reliably and was not used in my analysis. 

3.2.1 Beam-Beam Counter-BBC 

The BBC consists of two sets of quartz tube Cerenkov arrays, which 
measure relativistic charged particles produced in narrow cone around each 
beam axis (3.0 < r] < 3.9, 27r in 0). Positioned at 1.4 m from the PHENIX 
center point it has an outer radius of 30 cm and inner radius of 5 cm (see 
Fig. 3.2). 

Each BBC counter consists of 64 photomultiplier tubes (schematically 
shown in Fig. 3.3) equipped with quartz Cerenkov radiators in front. The 
dynamic range of each tube of the BBC allows to register 1-30 minimum 
ionizing particles which makes it the main Minimum Bias trigger detector in 
PHENIX for any collision species. 



Figure 3.2: Picture of BBC barrel be- 
fore the installation. 



Figure 3.3: Structural drawing of one 
BBC (beam view). Each box corre- 
sponds to one PMT. 



Each BBC PMT has intrinsic timing resolution of at = 50 ps and pro- 
vides high precision measurement of collision time and vertex position. For 
each collision, the BBC measures time of the collision with respect to the 
RHIC collider clock (synchronized with beam bunches). This time is usually 
referred to as the BBC t-zero and is one of the fundamental items of infor- 
mation about the collision in that it sets the start time for all the subsystems 
performing timing measurements. 

tQ^'" is calculated as a half sum of average hit time over individual BBC 
PMTs (denote as t^^*", tg^'^ for North and South detector). Obviously, vertex 
position can be calculated as a half difference of those variable. 

= {tr^ + tr^)/2; = (4^^ - tr')/^c (3.1) 

The vertex resolution in p + p collisions can be evaluated from the as- 
sumption that the multiplicity in the BBC is small and we have one hit in each 
BBC. Using Eq. 3.1, we can estimate in this case that az = crt/V^c ~ 1.2 cm. 
In central A collisions, accuracy of vertex measurement becomes much better 
{az < 0.3 cm) due to the averaging effect over all PMTs improves the timing 
resolution significantly. 

The BBC timing is an essential input to Level- 1 trigger providing online 
information about the vertex of the collision. Level-1 trigger electronics, gen- 
erate trigger accept signals if the vertex lies within ±50 cm from the PHENIX 
center point (to remove beam gas interactions and interactions inside the mag- 
net poles). This method enables the experiment to obtain a clean sample of 
events. Since the BBC trigger is efficient for a wide variety of interaction 
processes, it is referred to as a so-called "Minimum Bias" trigger. 
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3.2.2 Zero Degree Calorimeter 

The Zero Degree Calorimeter is a small area hadron calorimeter positioned 
^ 18 m from the interaction point along the beam axis. The main purpose 
of ZDC is measurement of the spectator neutron rate in nucleus-nucleus colli- 
sions. Spectator protons and other charged particles produced in the collision 
bend in the magnetic field of RHIC dipole magnets and miss the ZDC accep- 
tance. By measuring the flux of neutrons in heavy ion collision in correlation 
with the charged particles multiplicity in Beam-Beam counter, PHENIX ob- 
tains an information about the impact parameter on the collision (typically 
called the collision's "centrality" ) . A single ZDC counter consists of 3 modules 
each with a depth of two hadronic interaction lengths and read out by a sin- 
gle PMT. The ZDC provides timing and amplitude information similar to the 
BBC but with decidedly lower resolution. The energy resolution at the one 
neutron peak is approximately 21% [38, 39]. In Au+Au collisions the ZDC 
is an important part of Minimum Bias trigger, but in case of p + p collisions, 
presented in the analysis, it has no particular use since there are no spectator 
neutrons. 

Fig. 3.4 shows schematic view of ZDC detector. One can also see the 
projected trajectory for Au ions, spectator protons and neutrons at the ZDC 
placement plane A — A. 



m 50 
g 0, 

aJ-50, 

u 



DX 

dipole magnet ^ 



ZDCj 

7 A) 

I I I 1 1 I I I 

-20. -15. 







Intersection 
point 















-5. 0. 5. 

Meters 




ions 
neutrons- 



protons 



A-A 




-20. -10. 0. 10. 20. 30. 40. 

Centimeters 



Figure 3.4: Schematic view of ZDC detector. A) Top view of interaction 
region. B) Projected proton and neutron deflection area at ZDC placement 
plane. 



23 



3.3 Central Arm Detectors 

The Central Arm detectors provide unambiguous measurements of charged 
particle track information in mid-rapidity range {\ri\ < 0.35). The track's mo- 
mentum is determined by measuring its deflection angle in the magnetic field 
of the Central Magnet by the multi-wire Drift Chamber (DCH). A set 
of three Pad Chambers (PCI, PC2, PC3) helps to reconstruct the Z infor- 
mation of the track and remove the background tracks. The Ring Image 
Cherenkov detector (RICH) is essential for electron identification and pro- 
vides good e/n separation forpT < 4.8 GeV/c. Finally, the Electromagnetic 
Calorimeter (EMC) measures the energy, deposited by the charged particle 
or photon. The EMC also provides significant hadron/electron separation, 
crucial for electron analysis. 

The central arm subsystems are the most essential component of my elec- 
tron analysis and the final results rely on the stable performance of all of the 
described detectors. 

3.3.1 Central Magnet 

Although the Central Magnet (CM) is not a detector subsystem, central 
arm tracking relies on its stable operation. The Central Magnet consists of two 
coils (outer and inner) embedded into the massive steel yolk that generate an 
axially symmetric magnetic field in the region close to the interaction point. A 
picture of the CM during production is shown in Fig. 3.5 and the crossection 
schematic drawing is shown in Fig. 3.6. 

The main requirements to the Central Magnet fisted below: 

• Provide an smoothly varying magnetic field that can be mapped with 
0.2% precision. 

• Two coil operation, enabling creation of zero field at R — region in 
case of opposite coil polarization. 

• Minimal material in the PHENIX Central Area detector acceptance in 
order to minimize photon conversion background. 

• Magnetic Field strength should be significantly lower in the region of 
tracking detectors {R > 200 cm). This allows to assume a straight track 
model for tracking. 

• Movable configuration of the magnets in order to simplify detector as- 
sembling and commissioning. 




Figure 3.7: Magnetic field lines in PHENIX Magnet System. Drift Chamber 
location is shown by dashed box. 
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All the tasks listed were successfully accomplished, the magnetic field 
line contours in PHENIX Magnet system are shown in Fig. 3.7. Magnetic 
field strength as a function of radial distance R on Z = plane for different 
polarization of "outer" and "inner" magnetic coil is shown in Fig. 3.8. During 
Run02 PHENIX running period, used for this analysis, only the "outer" coil 
was energized which lead us to effective field integral for the charged track 
J Bdx = 0.78 [T ■ m]. An additional "inner" coil used in later runs helps to 
create an even stronger magnetic field in the acceptance in order to improve 
the momentum resolution for high pt tracks. 

The field reaches 0.096 T (0.048 T) at the DCH inner and outer radius 
at Z = and reach even smaller values at \Z\ = 80 cm which allows us to 
assume straight track in the DCH drift volume. 
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Figure 3.8: Magnetic field strength as function of radius for three different 
configurations of the magnetic coil polarizations. 
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3.3.2 Drift Chamber 

PHENIX Drift Chamber system contains a set of two "jet" -type multiwire 
detectors and is the main tracking device in PHENIX. It is placed at a radius 
of R E (202; 246) cm and consists of two identical arms each of which cover 
90°. The DCH performs the following tasks: 

• Accurate measurement of charged particle tracks in r-0 plane for deter- 
mination of their transverse momentum px- 

• Measure Z and 6 (inclination angle of the track with respect to Z axis) 
of the charged particle tracks together with PCI and BBC. 

• Provides input information for global tracking in PHENIX. 

During the construction of the Drift Chamber the following design re- 
quirements have been applied [37]: 

• Single wire resolution better than 150 fim in r-(f) direction 

• Single track reconstruction efficiency better than 99% 

• Two track resolution better than 1.5 mm. 

• Spacial resolution in the z direction better than 2 mm. 

In order to reach the specifications, the following design concept was im- 
plemented. A cylindrically shaped Ti frame was built as the support for wire 
nets with inner and outer radii of 202 and 246 cm and a length of 180 cm. A 
schematic view of one of the DCH arms is shown in Fig. 3.9. The active vol- 
ume of the DCH is filled with Argon(50%)/Ethanc(50%) gas mixture at STP. 
Charged particle tracking is done by a set of wire nets placed inside the gas 
volume of the chamber. There are total of 80 identical wire structures called 
nets around the azimuthal angle. Each net covers 1.125° in azimuth and is 
designed to measure the position of the track within its coverage by measuring 
the drift time {tdr) of the charge clusters, ionized by the incoming charged 
particle in the vicinity of the sense (anode) wire. 

The wire nets are subdivided into 6 separate sections along the radius, 
named in the following order (from inside to outside radius) XI, Ul, VI, X2, 
U2, V2. Each X wire net consists of 12 anode wires and measures the track 
trajectory in r — plane. Each U, V net has 4 sense wires and designed to 
reconstruct Z information for the track. Groups of 4 wire cells share the same 
electronics set and high voltage supply. Such a grouping is called a keystone. 
Fig. 3.10 shows the net configuration within one keystone. 
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Figure 3.10: 
Keystone. 



XI - plane 



Wire structure of DCH 



Figure 3.9: Schematic drawing of one 
Drift Chamber Arm. 

LA YOUT OF THE CELL 

Cell: New wire configuration 



I — I 213 

E 

u 

' — ' 212 
(0 

(0 211 



209 
20S 












f • 

— >* / 


ck g 

d 














d 
a 
d 


*U 


d 
b 
d 














a 

6 
e 


">■. 




o 

d 
O 

d 














<> 

d 

9 Fieic 

ft 


■I- It 


d 
d 

!>. 














o 
6 
d 

d 

6 

6 

d 




-I-H 
*« 

v-r'Anc 


o 
p 

<te d 








d 
o 

d 














o 
d 

6 




-Mt 


*\ 
d \ 














d 
d 
s 


'f, 

*» 


« 

g Ca 


:hod€ 












d 


■* 

« 


p 









X-axis [cm] 



Figure 3.11: Layout of the wire structure of XI celL 
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One cell consists of an anode wire net, surrounded by a pair of cathode 
nets. The region of ~ 2 cm between cathode and anode net is called drift 
region or charge collection region. Fig. 3.11 shows the wire structure of XI 
wire cell. 

The wire nets have a complicated wire and high voltage configuration in 
order to create a specially tailored electric field in the drift region. In total 5 
different values of high voltage are applied to different wires in order to create 
a narrow and isochronous alley in the drift region that supplies charge to the 
anode wire as shown in Fig. 3.11. The 5 voltages are named for the drift 
characteristic they control: 

• "Anode" (sense) wires - Read-out the charge, induced by drifting ion 
current. 

• "Cathode" wires - Create uniform electric field in the drift region 
(Vc « -4100 V). 

• "Back" wires - Block charge drift from one side of the anode wire to 
solve left-right ambiguity (Vb ~ —850 V). 

• "Gate" wires - Create a localized isochronous charge collection region 
and increase the field strength close to the anode wire {Vq ~ —2000 V). 

• "Field" wires - Separate individual anode wires drift regions and create 
a strong electric field around the anode wire {Vp ~ —2000 V). 

As a result, the electric field, created within one cell generates a well- 
localized charge collection region directed to one side of the net for odd an- 
ode wires and to the opposite side for the even wires. Fig. 3.12 illustrates 
GARFIELD [41] simulation of the regions which allow the ionized charge (dots) 
to drift to the anode wires. 

The drift velocity (vdr) within the drift region has a weak dependence 
on the electric field by the choice of gas mixture and typical field strength. The 
drift velocity as a function of electric field is shown in Fig. 3.13. The electric 
field is on the order of ~ 0.8 — 1 kV/cm in the area between the Gate and 
Cathode and is significantly higher between the Gate and Anode wire. The 
average drift velocity in the drift region is on the order of (vdr) ~ 50 /im/ns. 

The position of a hit within the cell can be calculated as using x — t 
relation: 

X = Vdr- {tdr - to) (3.2) 

where to is an important reference constant, corresponding to creation of 
the charge in the area of the anode wire. 
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Electron drift lines from a track 



Cell: New wire configuration Particle: 300 equally spaced points 

Gas: CjHg 50*. Ar 50*, T=300 K, p=1 atm 




x-Axis [cm] 



Figure 3.12: GARFIELD simulation of electric field lines inside XI cell of the 
DCH. Marked region display the charge collection zone of each anode wire. 
Circles represents charge clusters drifting towards the anode wire. 

Detection of the ionization signal would be impossible using room tem- 
perature electronics without gas amplification. Electrons travelling in the 
strong electric field can obtain enough energy between the collisions, to knock- 
out a secondary electron from a gas molecule, this secondary electron can then 
knock out another one and so on causing avalanche type multiplication of the 
charge. The threshold electric field is usually on the order of E^hr ~ 10 kV/cm 
and is reached very close to the anode wire. The electric field in the vicinity 
of the anode wire can be expressed as = ^ where Q is a charge per unit 
length in [1/] = Using GARFIELD we can simulate the distribution of 

the charges on each wire for given values of the potentials and wire geometry. 
From the test runs measurements (see Fig. 3. 14), it was found that Q = 290 V 
produce a single wire efficiency of 90 %. The standard voltage configuration 
was selected so that the charge on all planes exceed this limit. 

The DCH Front End electronics digitize the leading and trailing edge of 
the charge signal using the so-called "ASD-8" chip (Analog - Shaper - Discrim- 
inator) [40]. The thresholds for each input channel are set via an eternal DC 
voltage established by a DAG this is downloaded through an ARCNET net- 
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work. The typical running configuration uses a qthr = 6 fC threshold. The time 
of the leading edge is being digitized with a granularity of of the RHIC 
clock period or ^ 0.8 ns. , The trailing is digitized with twice coarser binning 
allowing the pulse width of each trigger to be measured. This measurement 
is used to reject especially narrow pulses as noise. The "Time Memory Chip" 
(TMC) functions by storing the continual running history of the leading and 
trailing edges for the previous 6 /xs in a circular memory buffer. Triggers force 
the chip to store one memory frame (an range from the past of user selected 
delay and depth) into one of 5 static memories. This allows the TMC chip to 
buffer hits from triggers that occur during the readout of previous data. This 
nearly eliminates deadtime from the PHENIX data collection system. The 
value of the "offset" (delay between the real data and the trigger arrival) is 
selected so that full drift time range could be read-out. Fig. 3.15 shows the 
typical time distribution shape from the DCH. The left edge of the timing 
distribution corresponds to the particles, depositing charge close to the an- 
ode wire^. The right edge of the timing distribution corresponds to the drift 
time of the electrons from the area of the Cathode wire. By measuring the 
half height time of the edges of the timing distribution we can make a rough 
estimation for and Vdr using Eq. 3.2. 



In the region between Gate and Anode wire we do not have Back side cancel- 
lation mechanism, this cause apparent double-counting of the signal and left-right 
ambiguity. 
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Figure 3.15: Time distribution for two DCH planes. Left and right edges fitted 
with error function [42]. 

DCH fine tuning 

It quickly became clear that in order to reach the design resolution, fine 
tuning of the calibration parameters (to and Vdr) need to be performed for each 
wire. This method was called internal alignment of the DCH and included: 

• Slewing correction - removal of to dependance of the width of the signal. 

• Vdr channel-by-channel alignment. 

• to channel-by-channel alignment. 

The most appropriate way of performing those type of corrections was 
based upon residual distributions. Residual by definition was denoted as 
At = to — [ti+t2)/2 where to is the time of the hit on the trial wire, ti and t2 are 
corresponding times for the neighboring wires hits. For convenience (as those 
calculations were performed on-line) the time unit for all time variables is going 
to be TMC time-bin (1 tb ^ 0.8 ns). Using the straight track assumption, it 
is clear that in optimal case At should be independent of all parameters and 
have a mean of zero. The side-standing wires (wires at the edge of the cell) 
have no neighbors and we can not calculate the residuals for them, they play 
a role of reference wires. 
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Figure 3.16: Residual distribution as a function of hit width before (left) and 
after (right) the slewing corrections. 




t„ Itb] 



Figure 3.17: Residual as a function of time fitted with linear function. 
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Figure 3.18: Comparison of Vdr for East 
XI wires obtained from the data resid- 
ual slope (circles) and GARFIELD sim- 
ulation (squares). 
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Figure 3.19: Distribution of the resid- 
ual distribution slopes (see Fig. 3.17) 
for all the DCH wires before and after 
drift velocity corrections. 
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Wc first look at the dependence of the residual as a function of the hit 
width (see Fig. 3.16). After including the fitted value of the mean shift as a 
correction to t-zero to = — dt'^''^^ {w) and performing at least three iterations 
we remove the dependence of residual on a hit width. 

The Drift velocity can be fine-adjusted using the field-off data. In this 
case, the majority of the tracks (not coming from decays or scatters) should 
go radially from the vertex with no defiection angle a ~ 0°. If we look at the 
distribution of the residual for particular wire as a function of time (shown 
in Fig. 3.17), the linear slope of the distribution that can only be caused 
by the difference of the drift velocities between the wire and its neighboring 
wires. It is possible to solve the corresponding system of linear equations 
for the corrections to the v^r variations with respect to the reference drift 
velocity. The solution require two constrains which can be selected arbitrarily. 
The best choice for choosing the constraints is by comparing the results of the 
corrections to the variation ofx — t relation slope from GARFIELD simulation. 
Fig. 3.18 shows the comparison of drift velocity profile obtained from data 
with the one produced from GARFIELD code. One can see nice agreement 
between expectation and the measurement. This method was performed on 
the plane-by-plane basics, the second order wire-by-wire corrections have been 
performed as a perturbation to the plane-by-plane corrections. Gomparison of 
the individual wire slopes after wire-by-wire fine tuning of the drift velocity 
is shown in Fig. 3.19. The final set of v^r correction coefficients was recorded 
as a multiplier for the global drift velocity, measured from the edges of the 
timing distribution. 

Fine adjustments to the calibration parameters were also performed on 
the to- Similar to the drift velocities, toS have strong dependence on the plane 
(probably due to the combined effect of the x — t relation non-linearity and 
the propagation time of the signal on on the ASD board). It also has a strong 
variation on the wire-by-wire basis due to possible geometrical displacement 
(driven by electrostatic sag) of individual wires and possibly different gas gain. 
The wire-by- wire variation of to was applied as an additional correction on top 
on plane-by- plane correction. 

In order to estimate the plane-by-plane variation of t^, the rising edge 
of the timing distribution was fitted at the constant level (1% - 5% of the 
maximum) as shown in Fig. 3.20. The results of the fit do not depend on the 
fit level (except for constant shift of the resulting time) and behavior is similar 
for both DGH arms (see Fig. 3.21). 

The local wire-by-wire to fine tuning was done plane by plane using the 
side-standing wires as a references and moving the middle wires to in order 
to zero the mean of the residual. The problem is easily reducible to a set of 




Figure 3.20: Determination of to from 
the slope of timing distribution at con- 
stant height for one DCH plane. Slope 
is fitted by Error function. 
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Figure 

3.21: Comparison of plane-by-plane to 
corrections for different DCH arms and 
different height cut. 



linear equations for dto j as a function of Ato «• 

After all the corrections are applied, the residuals distribution of the hit 
to the track reaches the design values of ~ 150/im which is shown in Fig. 3.22 
for both DCH arms. If we only look at the tracks far from the anode and 
cathode, this value can go down to ~ lOO/xm - which is probably the physical 
limitation due to the cluster arrival time smearing and the gas diffusion co- 
efficient. Needless to say that this value only indicates the relative accuracy 
of the tracking within one cell, data-based geometrical alignment of the nets 
within the DCH and matching to the outer detectors and vertex need to be 
done in order to perform global alignment of the Drift Chamber [42]. 
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Figure 3.22: Residual distribution of the hit to the track for East (left panel) 
and West (right panel) DCH arm fitted with double-gaussian [42]. 
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Tracking algorithm 

Tracking is the most critical step of the particle identification. Not only 
must the tracking operate with high single track reconstruction efficiency, but 
also it needs to be reliable in the high multiplicity environment of central 
Au + Au collisions and have a low "ghost" track rate. ^ The rate of the 
"ghost" tracks can be later reduced by association to other detectors (PC2, 
PC3, EMC). 

Much of the tuning to optomize the tracking performance was performed 
by me and I will summarize the key features of the method in this chapter. 

The ideal track (meaning 100 % single wire efficiency) should leave at 
least 6 hits in XI section and 6 hits in the X2 section of the DCH. In reality, 
inefficiencies cause the track to lose hits with a certain probability. The single 
wire efficiency can be calculated for each wire using the ratio of tracks, having 
a hit on this wire within some wide association window to the total number of 
tracks passing within that wire's active area. In order to remove the complica- 
tions of the cathode or anode crossing tracks, only the region confined close to 
the center of the cell is used for this calculation. The final efficiency map for 
p + p Run02 is shown in Fig. 3.23. One can see that the single wire efficiency 
is on the order of 95% for the middle wires and on the order of 90% for the 
side-standing wires (especially closest to mylar window of the chamber). This 
effect is well understood and was predicted from the Garfield simulations. It 
results from "edge effects" that reduce the field and hence the gain on the 
side-standing wires. 

Now we can estimate the tracking efficiency assuming that we consider a 
track to have at least N hits in XI and N hits in X2 for a given (constant) 
level of single wire efficiency. This can be exactly calculated using probability 
theory, the results of calculations for = 4, 5, 6 are shown in the Fig. 3.24. 
This plot show that having 95% single wire efficiency we already at the region 
of > 99% tracking efficiency if we chose N — A which was chosen as a minimum 
number of hits in XI or X2 layer for the track. 

The tracking in PHENIX is based on the assumption of the track having 
a straight line trajectory inside the Drift Chamber volume. First, the track is 
reconstructed in X ~ Y plane projection, determining a and (j) angles of the 
track, defined as indicated in Fig. 3.25. 

First stage of the trackfinding utihzes the "combinatorial hough trans- 
form" algorithm. Basic idea of the "combinatorial hough transform" is very 



^ Ghost tracks are pattern recognition solutions that did not actually come from 
real tracks. 
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Figure 3.23: Single wire efficiency map for p + p Run02. 



simple: we are looking for a straight tracks, having XI and X2 hits. Every 
straight line is described by two parameters (y=mx-|-b), both of which can 
be determined by any set of two points lying along the line. If we make all 
possible combinations between each XI and each X2 hit (laying within rea- 
sonable vicinity range from each other) we can calculate and histogram the 
line parameters. The line parameters (m,b) have infinite bounds and are im- 
practical to histogram, so we instead calculate, for each pair, the local angles 
Op and (f)p (same notation as on Fig. 3.25) to fill a 2-D histogram of ap vs. 
(pp. As the result - all the pairs of hits belonging to one track will create a 
localized peak on this histogram. The histogram can be replaced by 2-D array, 
which is called "hough array" and the local maxima in this array provide a 
"guess" value for a and of a track. In order to remove the possible bin-edge 
effects, a threshold is applied to the 3x3 bin sum around the local maximum 
of the "hough array" elements. This way we are able to clearly separate the 
background random combinations from the real tracks, to improve the initial 
guess parameters, we use weighted average of the neighboring array bins to 
improve the initial guess for a and (j). 




Figure 3.25: Single wire efficiency map for p + p Run02. 
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The threshold of the hough array is chosen Nf^r = 15, low enough to allow 
"4 hits in XI + 4 hits in X2" tracks to survive the cut. The binning of the 
hough array is also very critical issue: the coarser the binning, the more the 
probability to have all hits from one track localized to one bin. However if 
the bins are too coarse the accuracy for the determination of the initial guess 
parameters in decreased. Making the bins finer also has the drawback that 
hits from one track start to smear around the bins of "hough array" and the 
track may fall below the cut threshold. The bin size of the hough array was 
optimized by studying the reconstruction probability of one chosen track in 
high multiplicity Au + Au events (worst case scenario) . It does not matter at 
this stage whether we have a lot of background tracks, it is more important 
not to lose any. 

The next step of the trackfinding - is a "gentle" removal of the background: 

• Association of the hits to the track with (one hit ^ many track) corre- 
spondence. The association algorithm calculates the closest approach of 
the hit to the projected track guess and associates it if the hit is at least 
4 mm from the track guess. 

• Robust fitting of the track - iterative linear fitting using a weighting the 
hits, deweighting hits in accordance to their deviation from the mean of 
the previous iteration. This helps to remove randomly associated hits 
that are far off from the projected track. 5 fitting iterations performed, 
gradually rejecting mis-associated hits. Fig. 3.26 shows the track candi- 
dates (improved accuracy as compared to track guesses) after the first 
fitting stage. 

• The next step is sequential removal of the excess track candidates. 3 
stages of "gentle" removal perform the following procedure: 

1. Associate hits to the track with one hit ^ one track correspon- 
dence. In this step each hit will be associated only to the closest 
track. Association window is 4 mm for 1st and 2nd stage and 2 mm 
for the final association stage. 

2. Removal of the tracks that have less his than a threshold nthr as- 
sociated, nthr = 4 for the 1st stage, nthr = 6 for the 2nd stage 
and nthr = 8 for the final stage. This procedure gradually re- 
moves the tracks that have to few hits associated, returning their 
mis-associated hits to the real tracks. 

• Final fitting of the remaining tracks to the hits. Fig. 3.27 shows the 
tracks that were filtered to the output. Most all extraneous tracks are 
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removed by the method which proved to be extremely robust and efficient 
enabhng as to reach ^75% tracking efficiency in the most central full 
energy AuAu collisions. In low multiplicity environment ofp + p collision 
the tracking efficiency exceedss 98%. 




Figure 3.26: Event display snapshot with track candidates after the initial hit 
fitting. 




Figure 3.27: Event display snapshot with tracks on the final stage of track- 
finding. Same event as for Fig. 3.26 
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3.3.3 Pad Chambers 

The PHENIX Pad Chambers are unique background rejection devices 
using a novel pixelation scheme. Their primary purposes are following [37]: 

• Measurement of non-pro jective three dimensional spatial points, which 
are used for both momentum determination (p^) and pattern recognition. 

• Rejection of decays and photon conversion background at high px by 
tight matching requirements to the tracks measured by the DC. 

• Distinguishing electrons from other particles by accurate pointing of 
charged track to the RICH and EMC. 

• Charged particle veto in front of EMC. 

• Providing seed for tracks in charged high px Level-2 triggers and electron 
Level-2 triggers. 

The space points provided by PCS and EMC allow us to more accurately 
determine the track's actual trajectory through the RICH, an essential im- 
provement to the electron identification. All Z information for the track is 
obtained from PCI high precision Z measurement^ 

During the PC design of the following requirements were applied [44] 

• Very high efficiency (> 99%) and low occupancy(few % in most central 
Au + Au collisions) . 

• Good spatial resolution. 

• Low mass, in order to minimize secondary particle production and mul- 
tiple scattering. 

The PCs are multi-wire proportional chambers placed at radial positions 
of 2.5m, 4m, and 5m. Each detector contains a single plane of wires inside a 
gas volume bounded by two cathode planes. One cathode is segmented into 
an array of interlaced "pixel-pads". Each track fires three pixelpads. The 
coincidence reduces the false hit rate to be entirely negligible and localizes 
the track 3X better than a standard pixel chamber with the same number of 
channels. 

^The DCH can provide information about Z of the track using UV layers, but, 
the accuracy of this measurements is less precise than that provided by PCI alone. 
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A schematic view of the PC subsystem is shown in Fig. 3.28. The im- 
portant PC specifications achieved in RUN-2 are hsted in Table 3.2. The pad 
size for PCI is 0.84 cmx 0.845 cm to achieve less than 8% occupancy in most 
central Au + Au collisions. This gives a position resolution of 1.7 mm along z 
and 2.5 mm in r — (p. The pad size for PC2 and PC3 is chosen such that they 
have similar angular resolution compared to PCI. 




Table 3.2: Performance of Pad Chambers in RUN-2 [44]. 



Parameters 


PCI 


PC2 


PC3 


Pad Size (r-0 x z)[ 


cm^] 


0.84 X 0.845 


1.355 X 1.425 


1.6 X 1.67 


Single hit resolution 

{r-(f), z) in [mm] 


(2.5,1.7) 


(3.9,3.1) 


(4.6,3.6) 


Double hit resolution 

{r-(j), z) [cm] 


(2.9,2.4) 


(4.6,4.0) 


(5.3,5.0) 


Radiation Length 


[%] 


1.2 


2.4 


2.4 


Efficiency 


>99% 


>99% 


>99% 



42 



3.3.4 Ring Imaging Cerenkov Detectors 

The Ring Image Cherenkov (RICH) detector is the key component of 
PHENIX leptonic program. Not only does it have a nearly perfect rejection 
of pions over electrons up to ~ 5 GeV/c ( 1 x 10^^ error rate), but it also 
provides the Level- 1 electron trigger decision that enables us to collect rare 
electron and dielectron events. The main functions of RICH are: 

• Identification of electrons below px < 4.8 GeV/c. 

• Enable charged pion identification at > 4.8 GeV/c [42]. 

• Provide a fast Level- 1 trigger decision. In combination with EMC tile 
trigger, helps us significantly enrich the electron sample in high lu- 
minosity p + p collisions. Unfortunately, electron rate in Au collisions is 
too high to make the electron trigger effective. 

A schematic view of RICH detector is shown in Fig. 3.29. Each RICH 
detector has a volume of 40 m^. The spherical mirrors focus Cerenkov hght 

onto two arrays of photomultiplier tubes (PMT), each located on one side 
of the RICH entrance window. In order to achieve the design requirements, 
RICH performance has to satisfy the following specifications [37]: 

• c/tt separation at the 10^ level for single tracks. 

• The Photo Multiplier Tube (PMT) should have high single photon effi- 
ciency (> 99%). It should also have good timing resolution ( 300 ps) to 
reduce noise and contamination from albebo electrons generated during 
emc showers. 

• Minimal radiation length to reduce conversions inside RICH. 

The entrance of each PMT features a "Winston cone" of 50 mm diameter. 
The cone funnies light to the tube and increases the active area by reflecting 
light into the sensitive area of the tube that otherwise woTild have been missed. 
Each tube also has a magnetic shield that allows it to operate in magnetic field 
up to 100 Gauss. The radiator gas length seen by electron is 87 cm at = 
and 150 cm at 77 = 0.35, the average path length through radiator gas is 120 
cm. 

The RICH active volume is filled with CO2 gas at a pressure slightly above 
atmospheric. The gas has a Cerenkov threshold of '-/thr = 35 which is about 17 
MeV/c for an electron and 4.8 GeV/c for pion. The RICH can also be used 
to identify pions at pt > 4.8 GeV/c [42] Fig. 3.30 illustrates the principle of 
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Figure 3.29: A cutaway view of the RICH detector [45]. 
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Figure 3.30: Top view of the RICH and its optics. The tracks of the electrons 
and the emitted Cherenkov light cone are shown. Courtesy of Takashi Hachiya. 



Figure 3.31: Profile of tlie ring around tlie projected cliarged track intersection 
point. Courtesy of T. Kajihara. 



electron detection in the RICH. The Cerenkov photons generated by e~^,e~ 
and high momentum hadrons are reflected by spherical mirrors placed within 
the radiator volume. The photons are focused onto PMTs placed just behind 
the PHENIX central magnet. The pole tips of the magnet thus serve as hadron 
shields for the PMTs. 

Cerenkov photons are emitted at an angle of 6* = cos~^{l/{j3n)). These 
photons are focused as a ring of photons onto the PMT array, 

The total number of photo electrons for a charged particle above the 
Cerenkov threshold can be written as [9] 

Nnpe = L / ecedsin 6dE (3.3) 

TemeC^ J 

where ^"^^^a = 370cm~^el^~^ , L is path length of particles in the gas vol- 
ume, ec is the PMT Cerenkov light collecting efficiency and is the quantum 
efficiency of the PMT. 

Nnpe = Nlp^Lsin^O , where 

Kp. = (3.4) 

which quantifies the RICH electron detection performance. This number takes 
into account acceptance and quantum efficiency of the PMT and the property 
of the gas. In RUN-2, it is measured to be 116 cm~^ for CO2 gas. 

Each PMT has a diameter of about 2.5 cm, while the ring of photons 
reflected onto the PMT array has a radius of 11.8 cm. Fig. 3.31 shows the 
association are of the RICH ring with respect the the incident charged track. 
The RICH ring can be clearly seen as expected from the ring diameter for CO2 
gas. To reduce false hit rate, the number of PMTs for a given charged track 
are counted within 3. 4-12. 8cm from the projection - this value is called nO. 
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3.3.5 Electromagnetic Calorimeters 

The Electromagnetic Calorimeter (EMC) in the PHENIX is used to mea- 
sure the spatial position and energy of electrons and photons produced in 
Heavy Ion collisions. EMC also provides the means to trigger on rare events 
(high pt electrons and photons). The hadrons with kinetic energy more then 
200 MeV will not deposit significant energy in the calorimeter as the design 
and the thickness is dehberately chosen to be "uncompensated" . The detector 
covers the full Central Arm acceptance of —0.35 < r) < 0.35 and has two arms 
90° in azimuth each. 

The main specifications to the EMC design are listed below [37] : 

• Good energy and position resolution for electromagnetic showers. 

• Sub-nanosecond time resolution. 

• Comparatively low cost. 

To accomplish those goals, the basic design of EMC was selected to consist 
of 8 large sectors covering in total 60 of the PHENIX acceptance, 6 of the 
EMC sectors (E2,E3,W0-3) "Plumbum-Scintillator" modules (PbSc) and 
2 sectors of "Plumbum-Glass" modules (PbGl). The PbGl represents the 
greatest cost savings since the device was recycled from a previous experiment. 
The schematic view of the EMC sector is shown in Fig. 3.32. PbGl sectors 
of EMC were previously used by WA98 experiment and were re-installed in 
PHENIX. The PbSc sectors was built specifically for PHENIX in 1992 and 
designed as a set of lead-scintillator sandwich with readout by wavelength 
shifting (WLS) fibers penetrating the entire length of the detector cell (usually 
referred to as "tower"). The dimensions of one PbSc tower are 5.25x5.25x37.0 
cm^ and the effective depth of the EMC corresponds to 18 radiation length. 
The depth is chosen to obtain the optimal e/vr separation via E/p matching. 
For this analysis we use only PbSc sectors of EMC. Each tower measure the 
deposit energy of the electromagnetic shower in localized place on the detector 
surface. This enables us to look for the "clusters" of energy, locahzed to 
a particular block of towers and measure the total deposited energy under 
assumption of electromagnetic shower. The electron and photon leave all its 
energy in the EMC and localized to 2x2 towers with 85% probability, making 
the EMC an unique electron ID device in PHENIX. 




Figure 3.33: Design of PbSc EMC tower [37]. 
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The schematic design of one PbSc tower is shown in Fig. 3.33 there are to- 
tal of 18240 individual towers in PbSc EMC sectors grouped into 
25 "supermodules" each containing 64 towers. The module consists of a stack 
of alternating layers of 1.5 mm thick layers of lead, white reflective paper and 
4.0 mm thick polystyren-based scintillator tiles. Each stack has a drilled array 
of holes for the WLS fibers with 500 nm emission peak. The light is collected 
from fibers by a conventional PMT tube at the base of the tower. This design 
(so called "shish-kebab" type) prove to provide a perfect light collection and 
perfect time uniformity of the signal. The dynamic range of the PMT was 
chosen to perform the energy measurements starting from 0.1 up to 50 GeV 
with good linearity. 

The energy resolution of the PbSc EMC was measured on the test beam 
to be on the order of = 1.2 © ^-iid spatial resolution on the order 

of cr{x) = ^^1^ y-^ mm. Fig. 3.34 shows the energy deposited in the EMC by 

pion, proton and electron with different particle energy. One can clearly see 
very good separation of the electrons and hadrons (especially at > 1.0 GeV) 
which is the advantage of the electron identification capabilities of the EMC. 




Figure 3.34: Energy deposited in PbSc EMC by the pion, proton and electron 
of E = 0.5, 1.0 and 2.0 GeV. 
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Chapter 4 



Data Analysis 



This chapter describes the analysis procedure used to extract and iso- 
late the electron sample from heavy flavor semi-leptonic decays. Section 4.1 
describes the quality control, run and event selection used in the analysis. 
Electron identification cuts and their optimization is presented in Section 4.2. 
Section 4.3 explains the inclusive electron invariant crossection calculation. 
The estimation of the electron component from "photon" related decays of 
light mesons ("photonic" electron component) through the EXODUS Cocktail 
is summarized in Section 4.4. The final results for " non- photonic" electron 
crossection subtraction are presented in Section 4.5. In order to crosscheck the 
results of " non-photonic" electron component measurement, the independent 
"Converter subtraction" analysis was performed which is described in Sec- 
tion 4.6. Section 4.7 presents the Systematic Error estimations for inclusive 
electron crossection. Cocktail prediction and the subtracted "non-photonic" 
electron crossection. 



4.1 Quality Assurance and run selection 

The accurate run selection is essential for high precision measurement 
of Open Charm decay electron component. The contribution of "photonic" 
electron background is on the order of 80% at low Pt of total electron signal and 
even a small variation of total electron yield can cause a significant variation 
in the background-subtracted result. 

The other comphcation for the electron analysis is that we need to be 
certain that wc have uniform distribution of material in the acceptance. Any 
additional piece of equipment in the acceptance can cause a significant in- 
crease to creation of conversion electrons. Thus we need to apply an elaborate 
acceptance cuts in order to make a conversion rate uniform in the acceptance. 
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4.1.1 Acceptance cuts 

The acceptance for the electrons in PHENIX is best represented in terms 
of the track inchnation angle a and the azimuthal angle 0. The transverse mo- 
mentum of the particle is inversely proportional to a and can be approximated 
to the first order as pt ~ ^^^^ GeV/c. 

The Drift Chamber performance in the East arm was much more stable 
then that of the West arm and for this analysis we decided to use only the East 
arm acceptance. Applying very loose electron ID cuts (nO > 1, \dEMc\ < 5 ) 
we can plot the density of the electron candidates in a vs. space. Fig 4.1 
shows the electron acceptance of the East arm. One can see that big portion 
of acceptance is "shadowed" by conversions from Time Zero counter (TZR). 
This detector was installed into the PHENIX acceptance about 60 cm from 
the interaction point in order to improve the measurement of "start" time for 
the Time-of-Fhght detector. Unfortunately due to very large radiation length 
of TZR counter {Xtzr ~ 5.0%!!) it creates a very large rate of conversions 
far from vertex that creates a huge conversion background in the region of its 
shadow. ^ The region effected by the TZR counter is removed by the fiducial 
cut shown in the acceptance plot. Stripes on this figure depict acceptance 
holes for various PHENIX detectors. There is also a small portion of accep- 
tance (circled on the plot) affected by conversions from cable tray of the "New 
Trigger Counter" NTC detector, it is also removed by fiducial cut^. The small 
acceptance region at > 3.2 rad that is not affected by TZR "shadow" does 
not contain high momentum electrons which are of particular interest for Open 
Charm electron measurements. That is why analysis is using only EMC PbSc 
sectors E2, E3. 

Non-uniform conversion rate acceptance cuts are listed below: 

• 10 + 0.85 ■a\< 2.68; TZR counter cut 

• 10 0.85 • a\ > 1.95; NTC cable tray shadow 

• E2,E3; EMC sector cut 



Since the TZR detector debacle, any new detector placed in the PHENIX aper- 
ture has been required to submit a "detector impact statement" prior to its inclusion 

in our apparatus. 

^The NTC also failed to pass its "detector impact study" and is removed. Neither 
the NTC nor the TZR were ever used to produce a physics result in PHENIX. 



50 




Figure 4.1: Density of the electrons in a vs. space for East Arm. The 
area indicated by red arrows is the acceptance region used in the analysis. 
The various acceptance holes and additional photon conversion "shadows" 
indicated on the plot. 



The additional holes in the acceptance of the detector were studied start- 
ing from those closest to the interaction point [i.e. the DCH). Due to the 
bending of the track in the magnetic field, the azimuthal angle 0' at which 
track intercepts each detector component of the PHENIX will be shifted with 
respect to DCH 0, which is calculated at the "reference radius" -Rre/ = 220 
cm. The shift is proportional to a and is negative for interception with radius 
R > Rref and is positive otherwise, (pxi = (p + 0.06 ■ a corresponds to track 
angle in XI DCH plane, (f)x2 = 4' ^ 0.04 ■ a corresponds to track intersection 
of X2 DCH plane. By plotting track density in (pxi, 4'X2 coordinates we can 
clearly identify the dead DCH regions. 

The same analysis can be performed for the Pad Chamber (PCI) dead 
regions, (ppc can be approximated as (ppc = 4> ~ 0-13 • a. Fig. 4.2 shows the 
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track density in ^pc; -Z^ space. There is a PCI region of unstable performance 
that was cut-out. 
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Figure 4.2: Density of the electrons in Z vs. 0pc space. 



Final list of tracking fiducial cuts is summarized below: 



• Not{\^ + 0.06 ■ a - 2.562| < 0.005 & (Z < 0)) DCH dead region 

• Not(\<\) - 0.13 ■ a - 2.3651 < 0.025 & (Z < 0)) PCI dead region 

• Not{\<\) - 0.13 ■ a - 2.3251 < 0.025 & (Z < 0) & (Z > -15)) PCI dead 
region 



EMC dead area was calculated by photon density measurements in each 
EMC tower on a run- by-run basis [46, 47]. The final dead/noise map used in 
the analysis included all the towers that had a dead/noise flag set at least in 
one run. 3x3 tower region around the "bad" tower was fiducially removed in 
order to have more precise energy measurement in the vicinity of the "bad" 
tower. Fig. 4.3 shows the map of the dead towers in EMC that were removed 
from the analysis. 
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Figure 4.3: Dead/noisy EMC tower map for E2, E3 EMC sectors. 



4.1.2 Event selection 



The "Minimum Bias" event trigger in Run02 P+P was based on the coin- 
cidence of at least one hit each in the North and South Beam-Beam Counters 
(BBC). In order to keep up with high luminosity and to keep the constant 
bandwidth for Data Acquisition System, a trigger prescale logic was imple- 
mented. The live trigger rate was artificially reduced by only storing each one 
event out of Rscaie- Rscaie IS Called trigger prescale factor which depended on 
RHIC store luminosity and could be set to four possible values of 10, 20, 40, 
80. The Minimum Bias trigger was always prescaled. 
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In order to increase the rate of events containing energetic electron, 
PHENIX uses special Level- 1 electronics trigger called the ERT (EMC-RICH 
trigger). The basic principle of the trigger lies in the online summing of the 
energy signals in the EMC over a 2x2 tower region called a tile. If the signal 
from particular tile exceeds the tunable threshold value the specific bit is set 
in the data stream. To avoid edge effects, the tiles are overlapped and with 
2x2 summing the number of tiles equals the number of towers. The other 
bit is set once a RICH tile (4x5 tubes, overlapping) have a signal exceeding 
threshold. A spatial match between emc and rich tiles is an indication that 
the high momentum electron may have been detected in the particular region 
of detector.^ The trigger electronics issues the Local Level 1 (LLl) trigger 
decision for the PHENIX Global Level 1 system (GLl). The energy threshold 
of the ERT trigger can be adjusted by the threshold settings and was set to 
have a 50% registration probability for 800 MeV electron. The efficiency of 
the ERT trigger is discussed in more details in Section 4.3. The same elec- 
tronics can be used to fire on high energy photons by skipping the coincidence 
of the RICH bit. This trigger called "Gammal" was successfully used for 
high momentum 7r° measurements [47, 48] . The ERT trigger has a significant 
rejection power (the rate of the triggered events compared to Minimum Bias 
rate) Rert ~ 40 — 50, does not hmit the DAQ bandwidth, and requires no 
prescale factor. The proper normalization of the ERT trigger data should be 
done to the total number of live Minimum Bias triggers corresponding to the 
particular Run. It is crucial to not use the number of ERT trigger events 
for normalization as any noisy channel can artificially increase the rate of this 
trigger and a strong bias would be applied to the results. 

The collision vertex is measured by the Beam-Beam Counters (see Section 
3.2.1). The vertex resolution in p + p Run02 was 6z„t^ = 1.2 cm. Due to the 
specific geometry of PHENIX, the tracks originating from coUisions that are 
far from the center of the detector {Z^tx = cm) have a higher probability 
to interact with the material of the Central Magnet thus creating additional 
conversion electron background. The collision vertex distribution for the par- 
ticles that are primarily electrons (tight elD cut nO > 3) is shown in Fig. 4.4. 
One can see that the vertex distribution has an almost Gaussian shape with 
some additional structure for high Z^tx- In order to minimize the conversion 
background, we use tight vertex cuts for this analysis and only look at the 
events with \Zytx\ < 25 cm. 



^low momentum particles have displaced RICH and EMC tiles and fail the 
trigger. 
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Figure 4.4: BBC Z^t^ distribution for electron candidates (nO > 3). Bold line 
shows the vertex region used in the analysis. 

4.1.3 Run selection and event counting 

To filter the bad runs we looked at the distribution of all charged tracks 
{j)T > 0.4 GeV/c) with standard electron identification cuts (except for RICH 
nO > 1 cut, see Section 4.2) and all acceptance cuts. The ^ distribution for 
each run was normalized to the number of recorded Minimum Bias events 
Nmb- Then the ratio of the j^^^ distribution for given run to the same 
distribution for the chosen "reference" run (run having significant statistics 
and stable acceptance) was fitted with a constant R. Any significant deviation 
of the fit parameter R from one is an indication of additional dead area in the 
trial run. The criteria for the selection of the run was chosen to be: 

• R> 0.94 

• Xl < 2.0 

The run by run variation of R and xt presented on Fig. 4.5 (runs shown 
as red are considered to be "bad"). During Run02 we had an period with ad- 
ditional "photon converter" installed inside the PHENIX acceptance. "Con- 
verter" subtraction method explained in details in Section 4.6. Converter run 
period should be treated separately from Non-Converter run period. The total 
event statistics of the Converter and Non-Converter run periods summarized 
in Table 4.1. 
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Figure 4.5: Run-by-run variation of charged particle yield. 
Table 4.1: Statistics summary for Non-Converter and Converter run period. 





Nmb 


Nmb live 


Non-Converter "Total" 
Non-Converter "Bad" 
Non-Converter "Good" 


15 931 737 
540 061 
15 391 676 


475 849 920 
10 683 600 
465 166 320 


Converter "Total" 
Converter "Bad" 
Converter "Good" 


4 851 787 
361 847 
4 489 940 


264 284 240 
28 947 760 
235 336 480 
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4.2 Electron identification cuts 

Electron identification is one of tlie most critical parts of the analysis and 
required a precise tuning of the elD parameters. PHENIX is able to identify 
the electrons using the following parameters: 

• Number of RICH photomultipliers that have hit within the projected 
track ring - nO. 

• EMC matching - distance between track projection and EMC cluster in 
(p and Z coordinates. We denote (Iemc as (Iemc = \/ '^'^\mc + '^'^'emc- 

• Ratio of EMC deposited energy to particle momentum - E/p. 

Those variables have been used for all the current PHENIX electron re- 
sults [49, 50, 51]. All elD parameters were adjusted both for Data and Monte 
Carlo simulation to be identical and uniform (i.e. matching parameters have 
the mean value of zero and width of la independent of momentum and uniform 
throughout the detector acceptance). 

4.2.1 nO cut optimization 

The rejection power of separate elD cuts was tested by studying the E/p 
distribution of electron candidates before and after the cut. 

An initial assumption is chosen for the nO cut. This cut can not be 
set as low as one phototube since the random association background due to 
electronics noise would be too high. Thus nO > 1 was assumed to be lowest 
possible nO cut. The next step is to study what happens as the cut is tightened. 
Fig. 4.6 shows the effect of nO > 2 cut on the initial electron candidate sample 
(nO > 1) for pt > 0.4. The bottom left inlet shows the rejection power of 
the cut which shows the remaining portion of the particles after the cut. One 
can see a distinctive peak in E/p distribution which defines "real" electrons. 
The tail at E/p < 0.5 consists principally of random charged hadron tracks 
associated with real RICH hits and ofT-vcrtcx conversion electrons with mis- 
measured momentum. In order to estimate the effectiveness of the elD cut we 
check how much signal it removes in the E/p electron peak region and how 
strongly it cuts away the low E/p background. 

From the plot we can conclude that we lose ~ 10% of the electron signal 
and ~ 50 — 60% of the background. This loss of the electron efficiency is too 
large to afford and the rejection is not significantly high so we decided to use 
nO > 1 cut for the analysis. 
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Figure 4.6: a) E/p distribution for the electron candidates for different nO cuts 
b) Rejection power of elD cut nO > 2 in comparison with raO > 1 cut. 




Figure 4.7: E/p distribution for the electron candidates and rejection power 
for different cIemc cuts a) cIemc < 2cr b) (Iemc < 3cr c) (Iemc < 5cr. 



4.2.2 EMC matching cut optimization 

The similar studies have been done on the "adjusted" track matching 
dEMC = a/ d(j)%j^(j + dZ'^j^j(j parameter. dEMC variable cut has been tested 
for values of dEMC < 2, dEMC < 3, and dEMC < 5. The resulting rejection 
power of is shown in Fig. 4.7. 
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It is clear that (Iemc < 3 already cuts a big portion of low E/p elec- 
tron candidates leaving the peak statistics almost intact. The efficiency of a 
(Iemc < 2 cut at the peak is a; 80% which is a significant statistics loss. 
Therefore, (Iemc < 3 was chosen as an optimum for the analysis. 



4.2.3 E/p cut parametrization 

We expect electrons to generate an electromagnetic shower in the EMC 
and therefore register an energy equal to their momentum. The energy over 
momentum distribution not only enables us to identify electrons by also allows 
us to measure both the energy and momentum resolution. The resolution of 
E/p can be directly derived from and and can be written as: 



a{p) 



P 



= \/ <^lis + i<^DCH ■ py 



E V 



where is term due to the multiple scattering, {ancH ■ p) is DCH angular 
resolution, ac is a constant term of EMC energy resolution, - is an 

EMC energy resolution depending upon fluctuations in the number of particles 
produced in the EM shower. 

The fluctuations are independent and so 



^(-) ^ y + p4 "]I^C + <S+^ + i^DCH ■ P? (4.2) 

Eq. 4.2 was obtained assuming E ^ p ^ pt- From this equation one can see 
that at low Pt the main contributor to E/p resolution is energy resolution 
which is been overcome by momentum resolution term at high pt and starts 
grow linearly. 

The mean and sigma of E/p distribution for electron candidates was 
obtained as a function of pt by fltting each slice with Gaussian + expo- 
nential background (or Gaussian wherever background is negligible or non- 
exponential). Fig. 4.8 shows the flt results for different pt bins starting from 
0.4 GeV/c up to 5 GeV/c. 
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e/p , Pt > 2.0 && Pt < 2.5 | | e/p , p^ > 2.5 && p^ < 3.0 | | e/p , p^ > 4.0 && p^ < 5.0 




Figure 4.8: Fits to E/p distribution of electron candidates for different pt 
bins. 

Fig. 4.9 shows the mean and sigma distribution for the Gaussian compo- 
nent of the fit as a function of electron transverse momentum. One can see 
that we have a very clean electron sample with the background contribution 
slowly drifting to the lower E/p values going to higher momentum. 

From Fig. 4.8 one can see that our energy and momentum measurements 
are in good agreement. The apparent fall of the mean E/p is possibly related 
with the fact that at low momentum the inclination angle of the track becomes 
significant and EMC cluster starts to spread spatially and we start measure 
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Figure 4.9: a) Mean and b) Sigma of E/p Gaussian fit to electron signal as a 
function of pt (Sigma is fitted with Eq. 4.2 function). 



only a fraction of its total deposited energy^. 

The E/p resolution is fitted with Eq. 4.2 and the following energy and 
momentum resolution were derived: 

• ac® (Tms = (3.82 ± 0.86)% 

• (Temc = (8.47 ±0.28)% 

• crncH = (1.48 ±0.41)% 

This results are in good agreement with PHENIX measurements of energy 
and momentum resolution [47, 70] via other independent techniques. 

As the background level in p±p collisions is very low and significantly 
suppressed by EMC matching cuts, we decided to use a loose (±3cr) E/p cut 
for the analysis 

, E/p-<E/p> , o 
I aiE/p) I ^ 

Thus, the cut width follows the fitted sigma as a function of the track's mo- 
mentum. 

^The EMC cluster algorithm is tuned for photons that land at near-normal 
incidence. 
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4.3 Inclusive electron invariant crossection 

This section of the Thesis explains the procedure of single differential 
crossection calculation for single electrons starting from the "raw" distri- 
bution, correction of the "raw" electron distribution to full azimuthal & one 
unit in rapidity, estimation of background level, combining ERT and MB data 
sample, trigger bias correction of the final crossection and bin width related 
corrections. The final expression for MB and ERT inclusive crossection can 
be written the following way 

1 1 1 1 A/'e MB 1 O-pp tot ■ ^BBC 

Nmb 27r 2 pt Apr Ay erecoiPx) ^bias 

1 1 1 1 iVe ERT 1 CTpp tot ■ ^BBC 

Nmb live 27r 2 pr Apr Ay erecoiPr) ■ ^hias{PT) ■ (^ert{pt) 

(4.3) 



number of scaled minimum bias events in MB sample (Table 4.1) 
number of live minimum bias events in ERT sample (Table 4.1) 

" raw" electron count in pt bin 
rapidity range (±0.5 units in rapidity) 

reconstruction and acceptance efficiency (correction function) 
BBC trigger bias 

total p+p inelastic crossection [52] (42.2 ± 1.9) mb 
BBC efficiency for Minimum Bias [52] (0.516 ± 0.031) 
ERT trigger efficiency 



4.3.1 "Raw" electron yield 

A standard procedure for any spectroscopic measurement is to start with 
"raw" signal counting. First of all we selected an appropriate Pt binning which 
was chosen to match the bin boundaries of previous Au+Au single electron 
measurements [51]. Those bins that had significant statistics compared to 
Au+Au were split into two. The final choice for the binning is fisted below: 

• {0.4, 0.5, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6, 2.0, 2.5, 3.0, 4.0, 5.0} 



da 

^'Pt MB 
da 

^Pt ERT 



where 

Nmb 
Nmb live 

ApT 

Ay 

^recoiPT) 
^biasiPT) 
^pp tot 
^BBC 

^ert{pt) - 
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Figure 4.10: "Raw" electron count in pt bin for Minimum Bias and ERT 
trigger sample. 



Table 4.2: "Raw" electron count in pt bin for Minimum Bias and ERT trigger 
sample. 



Pt bin [GeV/c] 


Ne MB 


5Ne MB 


Ne ERT 


5Ne ERT 


0.4-0.5 


498 


22.31 


16 


4.00 


0.5-0.6 


278 


16.67 


69 


8.30 


0.6-0.8 


276 


16.61 


424 


20.59 


0.8-1.0 


104 


10.20 


846 


29.09 


1.0-1.2 


45 


6.71 


728 


26.98 


1.2-1.4 


19 


4.36 


433 


20.81 


1.4-1.6 


12 


3.46 


301 


17.35 


1.6-2.0 


10 


3.16 


263 


16.22 


2.0-2.5 


7 


2.65 


120 


10.95 


2.5-3.0 


4 


2.00 


30 


5.48 


3.0-4.0 








21 


4.58 


4.0-5.0 








5 


2.24 



Fig. 4.10 shows the electron statistics per bin for the ERT and Minimum 
Bias trigger data samples. Table 4.2 summarizes "raw" electron counting 
results. One can see that statistics in three highest px bins is quite low which 
is a limiting factor for Run02 single electron analysis. 
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4.3.2 ERT trigger efficiency 

The ERT trigger for Run02 was calculated for the J/ t/j analysis in p + p 
Run02 [49, 53] using single photons. The photon analysis uses only the EMC 
bit of the ERT trigger. However, the RICH trigger part the trigger does not 
introduce a significant efficiency loss. We can not use single electrons for 
the trigger efficiency measurement due to low statistics. The ERT efficiency 
calculation is trivial and described below. 

• Find a single photon cluster of energy E in EMC with tight identification 
cuts from a Minimum Bias event. 

• Check whether the ERT EMC trigger bit was set for this event and 
whether this particular photon fired the trigger. 

• The ratio of ERT registered yield dN,./dE to the total yield dNt/dE will 
give us ERT trigger efficiency eERT{E) 

Taking into account the fact that the electron momentum resolution at 
low pt is much better then energy resolution we use the trigger efficiency as a 
function of particle momentum instead of energy eEKrip)- 

The measured ERT trigger efficiency for the E2, E3 EMC sectors is shown 
in Fig. 4.11. The systematic error shows the maximum extent error of efficiency 
variation obtained by a 10% variation of the number of dead/noisy towers 
in the trigger simulation. The trigger efficiency may be underestimated for 
this analysis because we remove additional "bad" EMC towers (see Fig. 4.3) 
as compared to the analysis, but this difference is easily covered by the 
systematic error. Both the trigger efficiency and hi-lo limits of the systematic 
error were fitted with arbitrary functions, presented in Table 4.3 



Table 4.3: Fit results for ERT trigger efficiency and hi-lo systematic error band 
for E2,E3 EMC sectors. 





ERT efficiency 


Systematic error (hi limit) 


Systematic error (lo limit) 


E2 


1 


1 


1 


l+0.771/p4+o.473/p9 


l+0.451/p4+o.525/j)6 


l+1.307/;)4+o.685/p9 


E3 


0.96 


0.99 


0.94 


H-1.252/p4+o.588/pa 


l+0.681/p4+o.787/p6 


H-1.649/p4+l.435/p9 
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Figure 4.11: ERT trigger efficiency for E2, E3 EMC sector. The systematic 
error band to the efficiency from trigger simulation. 
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4.3.3 Hadronic background 

The major source of background for this analysis is a random coincidence 
of a charged track with a RICH cluster (thereby falsely identifying the track 
as an electron). This was studied previously in great detail for the Au + Au 
single electron measurements [50, 51] where this contribution is much more 
significant. In p + p collisions the multiplicity is low, thus the chance of ran- 
dom coincidence is significantly reduced (see Fig. 4.8 for background level 
estimation). The standard technique that is used in PHENIX offline soft- 
ware is so called "flip and slide" method which is based on creation of a fake 
("swapped") charged tracks by exchanging the North and South hits in all 
detectors except the drift chamber. This way we create an unbiased random 
track that then is being associated with outer PHENIX detector. The num- 
ber of RICH phototube that are associated with the " swapped" charged track 
denotes as snO. The distribution of E/p for the nO > 1 and snO > 1 is shown 
in Fig. 4.12 for Minimum Bias data sample. Unfortunately we can not use 
ERT trigger sample for those studies because the electron content of ERT 
events is strongly biased by trigger efficiency. One can see that the statistics 
of the purely random association of Minimum Bias sample is very small and 
alternative way to estimate the background contribution must to be found. 

We want to make an assumption at this point that the random association 
rate should not depend on the inclination angle of the track and, thus, it is 
not a function of the particles momentum. We would then find a constant 
probability erand thai the charged track is associated with a RICH ring. We 
can plot the E/p distribution for charged tracks and normalize it to the 
snO > 1 E/p distribution at low px- The normalization constant will be erand 
by construction. Fig. 4.13 shows the E/p distribution for nO > 1 , snO > 1 
, and charged hadrons scaled by erand = (3 ± 1.5{sys)) ■ 10~^ . In order to 
account for a qualitative comparison in the normalization, we put a large 
(50%) systematic error on this value. 

This probability does not include the effect of the ±3(7 E/p cut rejection. 
It adds an additional suppression of the random hadron association compo- 
nent by a factor 10 — 100. Fig. 4.14 shows the probability that randomly 
associated charged track "survives" the E/p cut. 



66 




0.2 0.4 0.6 0.8 1 1.2 1.4 



E/p 



Figure 4.12: E/p distribution for 
electron candidates (solid curve) and 
random association tracks snO > 1 
(dashed curve) for Minimum Bias 
events pt > 0.4 GeV/c. 
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Figure 4.13: E/p distribution for elec- 
tron candidates (solid curve), ran- 
dom association tracks (dashed curve) 
and charged hadron tracks scaled by 
^rand = 3 ■ 10""^ (thick solid curve) for 
Minimum Bias events px > 0.4 GeV/c. 




Figure 4.14: Rejection power of \E/p\ < Scr for charged hadrons. 
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4.3.4 (5-electron background 

5-electrons (also called "knock-on" electrons or (5-rays) refer to energetic 
electrons, that were knocked from the atomic shell of some atom in the de- 
tector volume. Depending on the construction of the detector and particle 
identification principles "knock-on" electrons may create a significant back- 
ground. 

In the case of PHENIX we need to take into account the rate of electrons 
emitted in RICH gas volume. The difference with previous effect is that delta 
electron, emitted with a reasonably small angle with respect to the initial 
hadron can create a hit in RICH and be misidentified as an electron. Thus, 
we need to estimate the probability for a hadron to emit delta electron that 
may fire RICH detector. To estimate the yield of 5-rays in RICH volume 
(100 cm of CO2) we used the formula ( 4.4, 4.5) [9]^ 

~ T^To (4-4) 



dTdx 2 A/32 T2 



where 



K - AnNArlmeC^ = 0.307075 [MeV cm^\ 

Z - atomic number of absorber 

A - atomic mass of absorber [gmol~^] 

z - charge of the incident particle 

T - kinetic energy of the electron 

7 - 7 the incident particle 

P - (5 the incident particle 

M - mass of the incident particle 

Tmax - maximal possible kinetic energy of the 5-electron (Eq. 4.5) 

RICH threshold for the electron is jthr — 35. Thus, minimal energy of 
electron that can "fire" RICH is = Ithr ■ itLsC^ ^ 17.9 MeV. Integrating 
Eq. 4.4 on T from Emm to T^ax we obtain the total yield of (5-electrons as the 
function of pion momentum (Fig. 4.15). 

The angle of the "knock-on" electron with respect to the incident pion 
can be calculated by the following formula [9] : 



^In calculations below we assume the validity of Rutherford crossection and 
spin-0 incident pion 
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Figure 4.15: Total 5-electron rate as a function of incident pion momentum. 



C0S(6') = {Te/Pe){Pmax/TTnax) (4.6) 

where 

Te,Pe - kinetic energy and momentum of the electron 

TmaxiPmax " maximum available kinetic energy and momentum of the electron 

Due to the RICH's geometrical acceptance, only when the 5-ray is pro- 
duced within cos{6max) > ^/n with respect to the pion direction will its 
Cerenkov radiation overlap in the RICH "ring" constructed around the pion 
projection point. The refraction index for CO2 gas n = 1 + 410 ■ 10~^ which 
gives the value for the maximum angle Omax = 28.63 mrad. From Eq. 4.6 one 
can calculate the minimal kinetic energy T^m of 5-electron which is deflected 
to an angle 6 = 6 max- Fig. 4.16 shows the range of kinetic energies for electrons 
that are emitted in ^ < Omax cone around the direction of incident pion. 

Now we can obtain the yield of ^-electrons that can be reconstructed in the 
RICH by integrating Eq. 4.4 from T^m to T^ax- (5-electrons rate as a function 
of incident pion momentum is shown in Fig. 4.17. One can see that we have a 
probability of es ~ 10^® for a pion to create such an electron (the rate is even 
smaller for incident particle of higher mass). The yield is significantly lower 
then previously calculated random association background of 3 ■ 10~^. Thus, 
(5-electron contribution can be neglected in PHENIX electron analysis. 
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Figure 4.16: The range of 5-electron kinetic energy that can be reconstructed 
by RICH, Train (sohd) and T^ax (dashed). 
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Figure 4.17: Rate of (5-electrons reconstruct ible by RICH as a function of 
incident pion momentum. 
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4.3.5 Acceptance Correction function 

In order to properly normalize the results of the measurements, we need 
to calculate the factor which describes the difference between an ideal An de- 
tector and PHENIX. This factor (see Eq. 4.3) is called the correction function 
(■recoipr) and takes into account the limited acceptance and track reconstruc- 
tion efficiency. 

The standard way to obtain the correction function is through the full 
simulation of particle detection probability assuming an ideally distributed 
input particle density. In the case of single electron analysis we " throw" single 
electrons, generated by the EXODUS [54] event generator, with the following 
input parameters 

• Uniform azimuthal angle distribution < < 27r 

• Uniform vertex Z^tx distribution^ \Zvtx\ < 25 cm 

• Uniform rapidity distribution —0.6 < y < 0.6 units 

• Uniform px distribution 0.0 < pr < 5.0 GeV/c 

The total statistics of our simulation sample was 3.98 • 10^ single particles 
(2.00 • 10^ positrons and 1.98 • 10^ electrons). 

The particles pass through the full detector simulation chain called PISA 
{PHENIX Integrated Simulation Application [55]). PISA is a GEANT-3 based 
simulation code that has been successfully used since 1992 to simulate real- 
istic particle propagation and detector response. The particle is "swimmed" 
through the tabulated Magnetic Field and GEANT-3 simulates the interac- 
tion of the primary particle with the material inside the PHENIX aperture. 
Both primary and secondary particles create a hit every time they enter the 
active area of the detector. This "Monte Carlo hit" information is stored in 
the "PISA output file". 

The next step in simulation is applying a realistic detector response to 
the MC hits. This procedure includes a smearing of hit position and timing 
with appropriate resolution, digitization of the timing and analog information, 
hit merging, applying time-of-fiight effects, and reproduction of registration 
efficiency & dead map for each detector subsystem, e.t.c. 

^ There is no strong dependence of reeonstruetion efficiency on Z^tx- This fact al- 
lows us to use a " flat" vertex distribution instead of realistic " Gaussian" distribution 
(shown in Fig. 4.4) 
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The final step of tlie simulation is the reconstruction of the simulated data 
using standard PHENIX offline code which is used for real data analysis. 

As an output we have a collection of reconstructed tracks. The offline 
software {evaluation package) maintains the relationship between MC track 
and reconstructed track. A "main contributor" scheme is used in that the MC 
track that provided most hits to any reconstructed track is then considered 
to be the "main contributor" to that track and is considered as the source 
of that reconstructed track. This helps us to obtain a direct correspondence 
between the reconstructed track parameters and the input track parameters. 

The first thing that needs to be checked in simulation is the momentum 
and energy resolution. Fig. 4.18 shows the mean and a of the difference be- 
tween reconstructed pt and initial pt mc transversal momentum as a function 
of Pt mc for electrons and positrons. One can see that there is a linear depen- 
dence of the S{pt) which needs to be removed from the simulation (we can not 
justify that this effect should exist in simulation and need to remove it and 
later treat it as a systematic error). 




Figure 4.18: Mean and sigma of difference between reconstructed and ideal pt 
in Simulation for electrons (left) and positrons (right). Linear fit is shown for 
mean distribution, sigma is fitted with Eq. 4.1 functional form. 
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The momentum resolution can be estimated from the fit to ct{pt — Pt mc) 
using functional form for the momentum resolution (Eq. 4.1). The Drift Cham- 
ber momentum resolution term in simulation is factor of two smaller then in 
the data: 



this means that we need to artificially worsen the momentum resolution 
in simulation and study what effect it may cause on the correction function 
(see Section 4.7). 

The momentum distribution of input particles is uniform and different 
from the ^ of real data. In order to take into account this difference, a 

weighting factor of {^^)Data dependent upon pr mc is applied to each Monte 
Carlo variable. This weighting procedure "artificially" adjusts the shape of 
the input MC momentum distribution to match the final Data momentum 
shape. The shape of the final inclusive electron distribution can be taken from 
Fig. 4.33 and we can assume the weighting function to be: 



The absolute scale is taken arbitrary in this formula as we are only interested 
in the shape of the input spectrum. Fig. 4.19 shows the comparison of the 
reconstructed tracks in Simulation weighted with w{pT mc) and reconstructed 
Minimum Bias (not final!) One can see that the shape of reconstructed tracks 
in the data agrees very well with Monte Carlo. 

Matching and elD cut parameters of the Simulations were adjusted in 
the same way as was done for the data. A comparison of the acceptance 
in the simulation and real data after applying the full elD cuts is shown in 
Fig. 4.20 for MB data sample (0.5 < pr < 2.0 GeV/c) and in Fig. 4.21 for 
ERT data sample scaled by ERT trigger efficiency (1.5 < pr < 5.0 GeV/c). 
The acceptance agrees well in all projections (0, (pEMC, Zemc)- 



• aMS = (0.87 ±0.05)% 



• (JDc/f = (0.74 ±0.02)% 




(4.7) 
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Figure 4.19: Comparison of distribution for weighted PISA simulation 



(circles) and MB data (squares). Ratio of in MB data to simulation 
(right). 
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Figure 4.20: Comparison of acceptance in 0, (pEMC, Z, Zemc for MB data 
(thin line) and weighted PISA simulation (thick line) for 0.5 < pt < 2.0 
GeV/c. 
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Figure 4.21: Comparison of acceptance in </>, (j)EMC} Zemc for ERT data 
(thin line) and weighted PISA simulation (thick line) for 1.5 < pt < 5.0 
GeV/c. 
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Now we have everything to calculate the correction function for simula- 
tion. By definition of the correction function: 



reco\PT Reco 



dN 

(iPT Reco 



■ w{p' 



'T Output) 



dN 

<^Pt Input 



w{pt 



(4.8) 



Input ) 



with 

d-VT Output 



where the ratio means the ratio of the histograms (for a given pt bin) filled 
{Pt Output) and ^^^^^ ■ w{pt input) correspondingly pt input 
denotes the transverse momentum of the input EXODUS particle, pt output 
is the transverse momentum of the "main contributor" PISA track associated 
to the reconstructed track with p^ Reco- This method of correction function 
calculation treats the weighting of the input and output distributions correctly. 

A rapidity cut of \y\ < 0.5 is applied to the input tracks in order to 
normalize the correction function to one unit of rapidity. 

The correction function for the e"*" + e~ simulation is shown on Fig. 4.22. 
The points are fitted in a range 0.4 < p^ < 4.5 GeV/c with a functional form 
that well-describes the shape. 
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Figure 4.22: Correction function erecoipx) for e"*" + e (full electron ID cuts). 
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The correction function indicates that we register 2.5% of the simulated 
particles in our acceptance almost independent on px. The apparent drop 
of erecoipr) at pt > 4.8 GeV/c is non-physical and is caused by the high 
momentum cut-off pr input < 5.0 GeV/c in the simulated particle sample. 

Correction functions for electrons and positrons separately are shown in 
Fig. 4.23 4.24. The shape of the correction functions for different charges is 
different because of highly asymmetric acceptance of the TZR counter " shadow" 
cut (see Fig. 4.1). 

x' I ndf 73.15/78 
pO 0.0240710.0007181 
pi -0.0132710.00178 
p2 -3.844 ± 0.2772 

p3 -0.9403 ±0.1 41 3 



"O 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

Pt [GeV/c] 

Figure 4.23: Correction function erecoiPr) ior e~ . Total correction function 
(dashed curve) shown for comparison. 

x' I ndf 61.14/78 
pO 0.023761 0.0008296 
pi 0.002495 ± 0.002656 
p2 -4.549 ± 0.7948 

p3 -1 .03 ± 0.3078 
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Figure 4.24: Correction function erecoipr) for e^. Total correction function 
(dashed curve) shown for comparison. 
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4.3.6 BBC Trigger Bias 

The remaining unknown in Eq. 4.3 is the BBC trigger bias ebiasiPr)- 
"BBC trigger bias" is PHENIX-specific term referring to the probabihty at 
which the BBC counter issues a Level 1 trigger decision for an event containing 
specific particle of interest. The overall BBC efficiency describes the fraction of 
the total p+p crossection registered by the BBC and was measured by Vernier 
scan to be tBBC = (0.516 ± 0.031) [52]. It is obvious that events with a hard 
parton scattering are more likely to be registered because the track multiplicity 
in the BBC is higher for these events. As an example, soft partonic scattering 
or worse still single- or double-diffractive scattering produce far fewer tracks 
in the BBC and are more likely to fail in generating a trigger. This means 
that of all events that contain a hard scattering process, the fraction recorded 
will be higher than the "inclusive" BBC trigger cross section. The fact that 
the trigger cross section depends upon the physics process is what we term 
"Bias". 

^biasiPr) was calculated [47] for 7r° production using the following tech- 
nique: 

• An unbiased sample of events was selected to be ERT 4x4 trigger with 
no BBC requirements. 

• 7r° was reconstructed through the 7r° ^ 7 + 7 decay in the EMC using 
the formula M^^ = 4 • • £"2 • sin'^{6^^/2). 

• The BBC trigger bias was calculated as a ratio of events with BBC vertex 
information reconstructed to the total number of ERT events. 

Fig. 4.25 shows BBC trigger bias as a function of neutral pion pt. One can 
see that the results agree with the value of e^ias — (0.75 ±0.02) independent of 
Pt [47] and, as expected, significantly higher than the inclusive BBC efficiency, 
esj^c = (0.516 ±0.031). 

This measured value of the constant BBC trigger bias is in good agreement 
with PYTHIA calculations of the BBC efficiency for hard pQCD partonic 
scattering processes [52]. Fig. 4.26 shows the PYTHIA simulation results 
for the BBC trigger efficiency as a function of the collision vertex, Z^tx: for 
different physical processes. One can see that the expected efficiency for pQCD 
hard processes is ~ 0.75 independent on the vertex position. Open Charm 
production should use this value of trigger bias since any collision process that 
can generate the charm quark's mass energy will certainly be a hard process. 
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Figure 4.25: BBC trigger bias for neutral pions as a function of 7r° px with 
the constant fit to the data [47]. 
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Figure 4.26: PYTHIA+PISA calculations for BBC efficiency as function of 
Z^tx for different physical processes. 



